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ABSTRACT 

A  semi-empirical  method  for  modelling  the  loss  of  electron  fluxes  in  the  earth's  mag- 
netosphere  was  developed.  An  equation  for  the  integral-energy  omnidirectional  electron  flux 
as  a  function  of  time  and  magnetic  field  strength  was  derived  from  pitch-angle  diffusion 
theory. 

This  flux  equation  was  the  basis  for  a  computer  data-fitting  program  written  at  the 
Air  Force  Weapons  Laboratory  (AFWL)  to  fit  the  AFWL  Trapped  Electron  Daoa  Base.  The 
program  utilised  a  least-squares  fit  and  incorporated  random  variations  of  the  characteristic 
exponential  loss  times  about  their  initial  values.  Ad  improved  table  of  initial  less  times  was 
compiled  for  use  with  the  program. 

The  derived  flux  model  showed  substantial  agreement  with  the  empirical  data  base. 
Representative  plots  of  computed  flux  over  raw  data  are  shown  for  L-values  of  2.3  and  2.4. 


VI 


I.  INTRODUCTION 


Objective 

Modern  technology,  both  civilian  and  military,  depends  greatly  on  the  use  of  satellite 
systems.  It  is  thus  important  to  the  Department  of  Defense  to  understand  and  predict  the 
environment  in  which  satellites  may  operate,  in  order  to  determine  survivability  and  vul- 
nerability requirements.  One  part  of  that  environment  is  the  electron  lluxes  which  may  exist 
in  the  trapping  regions  of  the  earth's  mngnetosphere.  Such  (luxes  may  arise  from  natural 
sources  or  from  injection  by  high-altitude  nuclear  bursts.  Because  such  fluxes  have  great 
capability  for  inflicting  damage  to  sensitive  satellite  components,  even  at  great  distances 
from  a  burst,  understanding  their  behavior  is  vitally  important. 

There  are  a  large  number  of  theoretical  treatments  in  the  literature  which  describe 
various  electron  loss  mechanisms  from  the  earth's  magnetosphere.  However,  the  empirical 
calculations  of  actual  electron  losses  are  often  conflicting,  and  frequently  are  restricted  in 
their  areas  of  coverage. 

Since  the  cessation  of  high  altitude  nuclear  testing  in  1962,  loss  studies  have  been  limited 
to  electrons  injected  into  the  magnetosphere  by  natural  magnetic  disturbances  such  as  solar 
storms.  Thus,  much  of  the  theory  which  has  been  developed  since  10R2  has  had  restricted 
opportunity  for  empirical  testing. 

The  objective  of  this  study  is  to  produce  an  improved  algorithm  to  calculate  losses 
of  electrons  from  the  magnetosphere.  The  algorithm  is  semi-empirical,  because  it  utilizes 
the  Air  Force  Weapons  Laboratory  (AFWL)  Trapped  Electron  Data  Base  (Ref.  37)  measure- 
ment of  electron  fluxes  (following  the  old  high  altitude  nuclear  bursts)  to  determine  the 
characteristic  loss  times  of  those  electrons. 

The  loss  times  thus  determined  and  the  algorithm  for  Qux  "decay"  will  then  be  incor- 
porated by  AFWL  into  the  SPECTER  computer  code  (Ref.  7-9)  for  improved  flux  calcula- 
tions of  satellite  environments. 


Scope 

This  study  is  limited  to  investigation  of  one  loss  mechanism,  pitch-angle  diffusion,  using 
a  theoretical  formalism  developed  principally  by  Shuk  (Ref.  47,48).  The  study  has  also  been 
necessarily  limited  to  those  regions  of  space  and  those  energies  covered  by  the  satellites  from 
which  the  data  base  was  complied.  The  loss  equations  thus  developed  in  this  study  are  only 
implicitly  functions  of  energy,  E,  and  L-coordinate  (explained  in  Chapter  II). 

The  sheer  volume  of  data  and  the  limited  time  for  this  study  also  limited  the  testing 
of  the  algorithm  developed  to  a  representative  subset  of  the  entire  data  base. 


Approach 

A  specific  solution  of  the  pitch-angle  diffusion  equation  is  derived  using  the  method 
proposed  by  Shulz  (Ref.  47).  This  solution  is  then  fitted  to  the  electron  flux  data  using 
a  least-squares  fitting  program  writteu  by  AFVVL,  with  the  coordination  of  this  author. 
The  solutions  are  performed  for  integral  fluxes  (above  threshold  energy  of  the  satellite 
detector)  and  for  constant  L-value.  The  solutions  are  thus  explicitly  functions  of  magnetic 
field  strength,  fi,  but  implicitly  functions  of  E  and  L  . 

Use  of  the  AFWL/NTCTS  computer  facilities  was  essential  to  thiB  study.  WH1  ovpr  72 
hours  of  actual  computer- processing  time  were  used,  ;vnd  over  1000  computer-produced  flux 
plots  were  examined. 


Presentation 

A  brief  Background  section  is  presented,  which  covers  3otne  physical  relations  and 
terminology  used  in  the  study  of  space  physics.  This  section  is  used  to  establish  uniform 
symbol  notation,  and  to  provide  a  common  starting  point  for  the  reader  unfamiliar  with  the 
specialized  language  of  space  physics.  It  may  be  skipped  by  the  reader  more  f;imiliar  with 
the  topic. 

Chapter  HI  outlines  the  theory  used  to  develop  the  flux  equations  which  were  fitted  to 
the  data. 

Following  the  theory  chapter  is  a  chapter  outlining  the  method  of  data  analysis. 
Descriptions  are  given  of  the  data  base  and  of  the  AFVVL  computer  codes  used  in  this  study. 

Representative  results  of  the  study  are  presented  in  Chapter  V,  in  both  tabular  and 
graphic  form.  Again,  the  sheer  volume  of  data  makes  inclusion  of  all  plots  and  lilted  values 
prohibitive. 

Conclusions  and  recommendations  for  improvements  for  further  study  are  presented  in 
Chapter  VI. 

The  Bibliography  includes  05  sources  and  provides  a  comprehensive  summary  of  existing 
literature  of  relevant  topics  about  radiation  trapped  in  the  magnetosphere  and  of  pitch-angle 
diffusion. 


H.  BACKGROUND 


The  equation  of  motion  of  a  charged  particle  in  magnetic,  electric,  and  external  gravita- 
tional fields  is: 

where 

r  =  r(t,r0,Uo);  -77  =  " 
m  =   particle  mass, 

g  =  gravitational  constant, 

<7  =  charge, 
H  =   magnetic  field, 
£  =  total  electric  field  (Rcf.  11:23-24.) 

In  a  uniform  magnetic  field,  a  charged  particle  will  move  in  helical  fashion  along  and 
around  the  field  lines  with  a  cyclotron  radius,  pe,  (also  called  gyroradius  or  Larraor  radius) 
about  its  "guiding  center".  The  particle's  motion  may  thus  be  separated  into  its  rotation 
around  the  field  line  and  the  motion  of  its  center  of  rotation,  or  guiding  center,  along  Hie 
field  line  (Ref.  11:24-25). 

The  period  of  cyclotron  rotation,  re,  is  defined  by  (Rcf.  39:5): 

2xm        2xpr 
qB  uj_ 

where  m  is  the  relativistic  mass: 


m  sss  m<(7  =  rru) 


Vi-KA2) 


and  where  Vj_  is  component  of  velocity  perpendicular  to  B.  The  cyclotron  frequency  is 
simply  2t  over  the  period  (Ref.  39:5): 

_  2x  _  qB 

tc  m 


The  angle  between  a  particle's  local  velocity,  v,  and  the  magnetic  field,  B,  is  its  pitch- 
angle,  a,  which  is  defined  as: 


(t)  =  *TCHv-f)- 


a  =  arccos 

In  a  uniform  static  magnetic  field,  vt  is  constant  and  v_[_  is  constant,  and  hence  a  is  constant. 
The  particle  will  then  move  with  a  uniform  circular  motion  around  the  field  line  and  uniform 
rectilinear  motion  along  the  field  line,  resulting  in  a  helical  motion  about  the  line,  ;.ince 

dp/dt     =     qv  X  B  =      ma, 

and 

mail       =     q(v  X  B),{         =     0, 


3!       =    0, 

yn        =     constant, 
and 

mS_L      =     ?(«  X  2?)j_, 
a  I         =     (qvj_B)/m       =     constant,  which  is  a  constant  centripetal  acceleration 

(Ref.39:4-7;  Ref.  11:23-38). 

The  concept  of  the  first  adiabatic  invariant  arises  naturally  from  the  gi  iding  center 
approximation.  The  first  adiabatic  invariant,  M,  (also  called  the  relativist  ic  magnetic 
moment)  is  deQncd  by: 

A 

\f  =  _  =  constant 

2m0B 

where  p_j_  is  the  particle's  perpendicular  momentum  in  the  guiding  center  ap  iroximation. 
The  assumptions  implicit  in  calling  M  invariant  are  that  the  spatial  variation  >f  B  is  small 
compared  to  pe  and  the  time  variation  of  B  is  small  compared  to  re  (Ref.  39:1  1-23). 

If  a  particle  has  a  constant  velocity  along  a  field  line  in  the  guiding  center  frame  of 

reference,  then 

■    '2    i    \  ■    2 

sin  a(s)        sin  a, 

~nTT  =  ~~n~~  =  con8tant 
B(s)  Bi 

where  s  is  the  arc  length  along  the  field  line  and  i  is  any  point  on  the  field  line.  This  assumes 
that  the  particle's  kinetic  energy  remains  constant  as  it  follows  the  field  line.  If  a  particle 
enters  the  field  at  point  /J,  with  pitch  angle  a,-,  its  velocity  along  the  field  line  is 


V||(a)  =  v  cosa(s)  =  v  •  \   1 -—  sin  a,-. 

If  the  field  is  increasing  in  the  direction  of  wj,  then  the  "mirror  point"  of  the  particle,  where 
its  parallel  velocity  is  zero,  is 

Bi      _      B(s) 

Mm.   =    7. —   7, 


sin'a,        sm"a(«) 
and  the  particle  has  a  local  pitch  angle  of  00°  at  that  mirror  point  (Ref.  39:34  42). 

When  the  magnetic  field  has  a  geometry  like  that  of  a  dipole  field,  increasing  3t  the 
poles  and  decreasing  in  the  midpoint,  then  the  particle  is  in  the  so-called  "magi  etic  bottle", 
and  is  trapped  between  the  mirror  points.  There  is  some  minimum  B  value  between  the 
mirror  points  which  is  called  Bo,  the  equatorial  field  strength.  Using  the  gi  iding  center 
approximation  and  approximating  the  earth's  magnetic  field  by  a  dipole  Gel  1,  a  particle 
trapped  in  the  earth's  field  has  three  distinct  motions.  It  circles  rapidly  about  :.  field  line,  it 
bounces  along  the  line  between  mirror  points,  and  it  drifts  slowly  in  longitude  around  the 
earth.  All  three  motions  take  place  with  different  speeds,  so  they  are  distinct  (Ref.  11:25; 
Ref.  29:34-65). 

From  the  mirror  point  definition,  and  from  the  velocity  equations,  we  see  that 

P*i 


M  : 

•2moB„ 


and  that 


W(*)  B(s) 


Note  that  the  particle's  mirror  point  depends  only  on  its  "injection"  point  and  its 
injection  pitch  angle,  not  on  its  initial  energy  or  velocity.  A  mirror  point  is  a  consequence 
of  the  field  alone,  and  all  particles  injected  at  the  same  point  "with  the  same  pit:h  angle  will 
mirror  at  the  same  point  on  the  field  line.  Of  course,  this  is  only  Lrue  if  no  external  forces 
are  acting. 

The  parallel  velocity,  v»,  of  the  particle  will  be  the  maximun  at  Bo,  wher?  B(s)  is  the 
minimum.  Thus,  a  trapped  particle  spends  most  of  its  time  near  the  mirror  po  nts,  and  the 
least  amount  of  its  time  transiting  the  equatorial  Geld  regions  (Ref.  29:34—14). 

The  bounce  period,  7&,  of  a  trapped  particle  is  generally  much  greater  than  its  cyclotron 
period,  rc,  and  is  defined  by: 

„   /"-■     ds  2  rm'  ds 

J.mi    «,(•)        vJ.mi    y/iZT{B{8)/Bm) 

•where  *mi  and  smi  are  the  mirror  points  on  the  field  line. 

As  a  particle  in  a  dipole  field  like  the  earth's  bounces  along  the  field  line,  it  also  drifts 
perpendicularly  to  the  field  line,  due  to  external  forces,  field  gradients,  and  fie'd  curvature, 
as  well  as  other  effects  such  as  time-dependent  field  changes.  This  drift  is  slew  compared 
to  the  bounce  period.  As  the  particle  bounces  and  drifts,  it  traces  out  a  surface  between  its 
mirror  points  and  around  the  earth,  called  a  drift  shell  (Ref.  39:9-19). 

The  concepts  of  the  second  and  third  adiabatic  invariants  arise  from  the  above  behavior. 
If  the  forces  acting  on  the  particle  remain  almost  constant  over  its  bounce  period,  the  second 
adiabatic  invariant,  J ,  is  defined  by 


=  j  p„  ds  =  4  /       PQ  ds 


where  pj|  is  the  momentum  component  parallel  to  B  and  ds  measures  arc  len^h  from  the 
equator  (Ref.  45:12;  Ref.  6:3-31).  If  the  forces  acting  on  a  particle  remain  almost  constant 
over  its  drift  period,  the  third  adiabatic  invariant,  <J>,  which  is  the  magnetic  flux  enclosed  by 
a  drift  shell,  is  defined  by 


4>=  f 


Aq  *  ds 


where  Aq  is  the  magnetic  vector  potential  and  where  the  integration  is  performed  over  a 
curve,  s,  which  lies  in  the  drift  shell.  The  third  invariant,  <1>,  is  defined  and  c  imputed  for 
the  drift  shell  of  the  guiding  center,  with  constant  field,  and  not  for  a  drift  ?hell  which  a 
particle  may  physically  trace  out  under  short-term  conditions  (Ref.  45:12;  Ref  39:76-79). 

The  guiding  drift  shell  of  a  particle  may  also  be  referred  to  as  an  "invariant  surface," 
which  is  composed  of  field  lines  which  end  at  the  mirror  points.  The  three  adiaba.  ic  invariants 
uniquely  define  an  invariant  surface  (Ref.  6:3-31  to  3-32). 

The  more  common  set  of  parameters  used  to  define  a  particle's  position  in  the  mag- 
netosphere  (or  to  define  an  invariant  surface)  is  the  B-L  coordinate  system.  The  8  parameter 
is  magnetic  field  strength,  and  L  is  the  Mcllwain  L  parameter  defined  by  L  =  iq/  He  where 
RE  is  the  earth's  radius  (?a  6371  km)  (Ref.  8:3-33)  and  (Ref.  39:53). 

The  value  of  L,  therefore,  is  equal  to  the  distance,  in  earth  radii,  of  t!  e  equatorial 
point  on  a  Geld  line  in  a  dipole  field.  If  the  field  is  not  symmetric  (not  a  pe:  feet  dipole), 
the  invariant  surface  is  not  so  well  defined.  However,  for  most  field  lines  iu    .he  trapping 


region  of  the  earth's  magnetosphere,  L  varies  by  less  than  one  percent  along  the  line,  so 
the  D-L  system  is  adequate.  The  less-than-exact  symmetry  of  the  earth's  dip  Die  field  and 
outside  forces  may  cause  particle  drift  to  deviate  from  a  perfect  azimuthal  -ourse.  This 
means  that  its  L  value  will  vary  in  the  course  of  a  drift  period.  This  variatioi  of  L  values 
is  not  significant  below  L  ==  3,  but  may  be  so  above  that  value.  The  average  //-value  of  all 
intersecting  invariant  surfaces  is  called  an  L-sheli  This  L- value  defines  a  set  of  surfaces  along 
field  lines,  which  may  end  at  dillering  mirror  points.  The  //-shells  are  considered  to  intersect 
the  earth's  surface,  even  though  mirror  points  do  not  extend  through  the  atmo:  phere  to  the 
surface  (Ref.  6:3-33  to  3-35;  lief.  39:53). 

The  various  mechanisms  which  operate  to  cause  trapped  particles  to  change  /./-shell  or 
to  be  lost  from  the  trapping  region  will  necessarily  violate  one  or  more  of  the  adiabatic 
invariants.  For  example,  pitch-angle  diffusion  violates  M  or  /,  or  both;  and  radial  diffusion 
violates  <b  (Ref.  45:48).  These  concepts  are  discussed  further  in  Chapter  111. 


EI.  THEORY 

While  the  motion  of  particles  trapped  in  a  magnetic  field  is  well  understood,  the 
mechanisms  of  losses  from  magnetic  fields  such  as  the  earth's  magnetosphere  are  less  well 
understood.  It  is  generally  agreed,  however,  that  pitch-angle  diffusion  into  the  loss  cone  is 
one  of  the  predominant  mechanisms  for  removal  of  charged  particles  for  mid-range  /^-values. 
The  loss  cone  angle  is  the  lower  limit  of  pitch  angle  for  trapped  particles.  Any  particles  with 
smaller  pitch  angles  will  mirror  in  the  sensible  atmosphere  and  will  be  lost  by  atmospheric 
scattering. 

Roberts  (Ref.  41)  has  noted  that  a  pitch-angle  scattering  mechanism  mest  be  extant 
for  pitch-angle  diffusion  to  occur.  Such  a  mechanism  would  necessarily  violate  one  or  more 
of  the  adiabatic  invariants. 

It  is  not  the  purpose  of  this  study  to  ascertain  the  true  source  mechanism  of  pitch- 
angle  diffusion.  Such  a  study  is  beyond  the  scope  of  this  paper.  Indeed,  the  search  for 
this  mechanism  has  been  going  on  for  decades.  Roberts  (Ref.  41)  has  postulated  that  the 
"perturbation  forces"  causing  such  diffusion  may  result  from  turbulent  ambient  electric  or 
magnetic  fields  or  from  collisions  with  other  trapped  or  non-trapped  particles.  Lyons  (Refs. 
20-28)  has  formulated  extensive  theory  based  on  resonant  interactions  of  so-calli  d  "whistler" 
VLF  waves  with  the  trapped-particles'  gyrofrequencies:  "cyclotron-resonance".  Regardless 
of  the  actual  physical  mechanism  behind  pitch-angle  diffusion,  its  general  treatment  is 
mathematically  the  same,  and  the  physical  results  are  the  same.  (Ref.  41:308)  This  study 
follows  the  methods  of  Roberts  (Ref.  41),  Shulz  and  Lanzerotti  (Ref.  45)  and  Shulz  (Ref. 
47)  in  assuming  that  pitch-angle  diffusion  is  an  operative  process.  One  of  the  purposes  of 
this  study  is  the  formulation  of  the  equations  necessary  to  validate  that  assumption  against 
experimental  data,  and  to  perform  that  experimental  validation. 

It  should  be  emphasized  that,  while  pitch-angle  diffusion  is  assumed  in  this  paper  to 
be  the  dominant  loss  mechanism  for  L-values  and  altitudes  considered,  it  is  not  the  only 
loss  mechanism.  Several  investigators  have  formulated  radial  diffusion  (cros:;-L)  theories 
(for  example,  Walt  (Ref.  58)  or  Tomassian  (Ref.  54)),  coulombic  or  collisioi-scattering- 
diffusion  theories  (for  example,  Wentworth  (Ref.  60))  and  multiple  diffusion  theories,  such  as 
combinations  of  pitch-angle,  energy,  and/or  radial  diffusion  (for  example,  Walt  (Ref.  58)  and 
Lyons  (Ref.  26)),  all  of  which  show  some  agreement  with  experimental  data.  la  particular, 
atmospheric  scattering  is  obviously  a  dominant  force  at  very  low  L  values  (Ref.  30). 

The  introduction  of  multiple  loss-mechanisms  makes  explicit  solutions  of  iny  diffusion 
equation  extremely  difficult.  To  simplify  the  problem  and  render  it  amenable  to  the  method 
of  Shulz  (Ref.  47),  this  study  assumes  that  "as  a  rule,  radial  diffusion  enables  the  radiation 
belts  to  become  populated  from  an  external  source  (or  rearranges  particles  injected  by  an 
internal  source),  while  pitch  angle  diffusion  causes  particle  loss  to  an  atmospheric  sink"  (Ref. 
45:48). 

Pitch-angle  diffusion,  while  simple  in  concept,  is  complicated  in  detail.  Roberts'  (Ref. 
41:307-337)  treatment  of  the  general  mechanism  is  particularly  descriptive: 

At  the  magnetic  equator,  a  particle's  pitch  angle,  ao,  is  determined  by 
(1)  x  =  COSfto  = 


V<  v + pj-.!) 


where  F\0  and  Pj_„  are  the  parallel  and  perpendicluar  components  of  the  particle's 
momentum  with  respect  to  the  Qeld  line  at  the  equatorial  value  of  Bq.  Roberts  next  defines 


(2) 


Wi_  = 


:m0 


"i.-* 


2  TTlfl 


where  mo  is  the  rest  mass.  These  are  defined  since  the  first  adiabatic  invariant,  M,  is 
proportional  to  Wj_0: 

(3) 


and,  if  Wj|„  <  Wj_0,  the  second  adiabatic  invariant,  /,  is  approximately  pre portional  to 

W  J  ~  "T" 


where 


rg  =  the  particle's  bounce  period, 
relativistic  mass  m 


rest  mass 


mo 


A  particle's  path  may  be  defined  by  plotting  its  values  of  Wj_0  and  Wj|0  as  in  Figure  1. 
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Fig.  1.    Loss  of  a  trapped  particle  by  a  random  walk  into  the  pitch- 

ande  loss  cone.  (Ref.  41:308) 

Note  that  lines  of  constant  equatorial  pitch  angle  pass  through  the  origin,  because 


tan   ao  = 


-Wj. 


w». 


Any  force  which  violates  either  M  or  /,  the  first  or  second  adiabatic  invariant,  will 
produce  "diffusion"  of  the  particle's  pitch  angle  as  shown  in  Figure  1  as  Wj_0  and  W^0 
change  with  time.  When  a  particle's  equatorial  Ditch  angle  diffuses  to  the  vali  e  of  the  loss 
cone  angle, 


(5) 


x.  —  cos  a,,  =  \/l  - 


Jk 


00 


where  Bioo  is  the  value  at  100  km  altitude,  then  it  is  lost  by  atmospheric  scattering.  "The 
loss  cone  serves  as  a  'sink'  for  particles"  undergoing  pitch  angle  diffusion  (Rcf.  41:307). 

Roberts  makes  two  very  important  points  about  pitch-angle  diffusion  as  a  oss-mechan- 
ism: 

(1)  with  no  source,  the  entire  radiation-belt  would  be  depleted  of  particlcr, 

(2)  the  loss  cone  approximation  is  just  that,  since  the  atmosphere  is  not  sh:.rply  defined 
at  100  km.  However,  if  atmotpheric  pitch-angle  scattering  is  not  the  primary  focus, 
the  loss  come  approximation  may  be  useful.  "Naturally,  when  the  loss-  one  concept 
is  used,  detailed  agreement  between  theory  and  experiment  cannot  be  expected  in 
the  region  near  the  edge  of  the  loss  cone"  (Ref.  41:308). 

Both  Roberts  (Ref.  41)  and  Shulz  and  Lanzerotti  (Ref.  45)  develop  the  pitch  angle 
diffusion  equation  from  the  Fokker-Planck  equation  using  a  particle  distribution  function, 
/• 

The  Fokker-Planck  equation  is  a  formalism  which  arose  from  the  study  of  Brownian 
motion,  and  which  is  used  frequently  in  transport  theory  (Ref.  42:308).  The  characteristic 
Fokker-Planck  equation  for  trapped  particles  (which  ignores  radial  diffusion)  i: 


dJ. 

dt 
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x  T[y)  dx 


T{y)D»  y~ 
ox 
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where  the  first  term,  subscripted  v,  represents  non-stochastic  (mean)  energy  loss  to  the 
atmosphere,  the  second  term  represents  pitch-angle  diffusion,  and  the  third  ter  n  represents 
range-straggling  (energy  diffusion).  The  term  T{y)  is  defined  below,  and  the  D^E  aQd  Dxs 
terms  are  the  characteristic  diffusion  coefficients;  f  is  relativistic  mass  ratio  (R<  f.  45:55-58). 

Roberts  uses  a  distribution  of  particles  in  a  "tube"  of  force  about  a  fi<  Id  line,  and 
Shulz  and  Lanzerotti  use  a  phtue-tpace  density  dutribution  function  which  s  essentially 
equivalent  to  Roberts'.  Since  this  study  follows  the  methods  of  Shulz  (Ref.  4  ')  and  Shulz 
and  Lanzerotti,  a  discussion  of  phase  space  is  necessary  (Ref.  45:15-22). 

Any  moving  particle  can  be  described  by  specifying  its  three  position  coo-dinates  and 
its  three  canonical  momentum  components.  This  completes  a  six-dimension  "i>hase-spacc" 
in  which  a  particle  can  travel  in  time. 

If  there  exists  a  system  of  a  large  number  of  particles  in  phase-space,  th<  system  can 
be  described  by  a  six-dimensional  distribution  function  j{Pi,qi,t)  where 


Pi  (t  =  1,2,3)    are  canonical  momentum  components, 
7,  (t  =  1,2,3)    are  position  coordinates,  and  (  is  time. 


Thus,  /  d^Pd^Tf  is  the  number  of  particles  instantaneously  occupying  6-D  volume 
d3Pd3q.  According  to  Liouville's  theorem,  the  phase-space  volume  containing  the  system 
of  particles  moves  incompressibly  through  phase-space. 

Since  f  is  an  awkward  quantity  to  deal  with  physically,  Shulz  and  Lanzerctti  note  that 
P  transforms  to  the  more  familiar  p  as 


(8) 


>-»+(!)* 


where  A  is  the  electromagnetic  vector  potential. 

Hence,  f{?,q,t)  =  f(p,?,t)  since  q  =  ?  and  the  P  to  p  transformation  has  a  unit 
Jacobian.  The  position-momentum  distribution  function  f[p,r,t)  defines  the  particles  oc- 
cupying the  6-D  volume  d3pd3r  at  any  given  t.  The  pitch  angle  diffusion  equation  is  given  in 
terms  of  a  phase-space  distribution  function  which  is  numerically  equal  to  the  more  easily- 
definable  position-momentum  distribution  function  (Ref.  45:15-22).  At  constant  energy  and 
L-shcll,  and  under  the  action  of  some  source  S,  the  diffusion  equation  in  phase- >pace  can  be 
written 


{') 


dt 
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x  T(y)  dx 


*T{y)D„ 


dx 


where  T{y)  is  the  quarter- bounce  integral  path  length  function  of  y  =  y/T  -  x2  where 
x  =  cosao-  The  quarter  bounce  integral  gives  the  length  of  the  trajectory  of  ;i  particle  (in 
units  of  Ro,  the  distance  to  the  equatorial  crossing)  from  the  equator  to  the  mirror  point. 
The  exact  definition  of  T  is 

T.r— 

Jo       cos  a 
where  n  is  the  distance  along  the  field  line  and  a  is  the  local  pitch  angle  (Ref.  10:4029—1030). 

Dxs  is  the  bounce-averaged  diffusion  coefficient  as  a  function  of  i  (Ref.  41  and  Ref. 
45).  The  boundary  conditions  are  that  /  =  0  at  x  =  xe  (some  cutoff  value)  and  f{x  =  0)  is 
finite. 


Shulz  and  Lanzerotti  (Ref.  45:182)  rewrite  the  above  as 
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and  further  approximate  that  the  second  term  is  negligible  for  x  *Z  1  since  i2  +  y2  =  1. 
If  T[y)  «  T{\),  the  second  term  disappears  and  the  equation  is  a  diffusion  equation  in 
cylindrical  coordinates.  (These  approximations  are  not  used  in  this  work.  This  development 
is  used  to  aid  in  understanding  the  exact  solutions  which  follow  from  Shulz  Ref.  47)).  If 
Dxs  and  5  are  then  independent  of  x,  and  if  f(x,  t)  =  X(x)  T(t),  the  eigenfurctions  of  (8) 
are  Bessel  functions  of  order  zero.  The  general  solution  to  (8)  would  then  be 


(«) 


7(M)  =  /«,(*)  +  E  a«W  J°  (*»  ~l) 


where  /ro  's  the  steady-state  solution 


(10) 
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The  Kn  are  the  teroes  of  /<>  (n  =  1,2,...)  and  the  an(t)  vary  as  e~''/r"',  where  r„  = 
x*/(D2k%).  Thus,  the  pitch  aagle  distribution  (and  hence  the  directional  or  omnidirectional 
flux)  is  shown  to  be  the  sum  of  a  steady  state  and  higher  order  eigenmodes  (Ref.  45:160-168). 

The  steady  state  can  be  thought  of  as  the  "normal"  or  "quiet-time"  value  of  flux  which 
exists  in  equilibrium  with  the  source  5.  An  "injection"  of  particles  would  theo  result  in  a 
perturbed  distribution  function  with  several  eigenmodes,  each  with  a  characteristic  decay- 
time  r„.  The  higher  order  modes  decay  faster,  and  eventually  only  the  fundamental  mode 
would  remain,  which  would  decay  exponentially  to  reach  the  steady  state  (after  "inCnite" 
time)  (Ref.  41  and  45).  There  is  ample  evidence  that  such  a  process  does  indeed  occur. 
Roberts  (Ref.  41)  cites  Explorer  XV  data  for  the  28  October  1962  Russian  explosion  as  an 
example.  Rosen  and  Sanders  (Ref.  43)  also  note  that  decay  is  faster  immediately  after  solar 
magnetic  storm  activity  than  during  quieter  periods.  The  similarities  in  Figures  2  and  3 
illustrate  the  way  in  which  a  decay  of  higher  eigenmodes  (Figure  2)  can  approximate  the 
temporal  evaluation  of  an  actual  electron  distribution  (Figure  3). 

The  primary  difficulties  with  the  above  discussion  arise  from  the  approximations  that: 

(1)  T(y)  =  T{\)  =  constant 

(2)  x2  <£  1  (orywl) 

(3)  Dxs  varies  little  (or  not  at  all)  with  x. 

Roberts  (Ref.  41:31 1)  uses  the  argument  that  the  full  spiral  path  (and  hence  T[y))  varies 
only  by  a  factor  of  1.4  as  x  varies  from  0  to  0.9  and  by  1.9  as  x  varies  from  0  to  1.  Shulz 
and  Lanzerotti  (Ref.  45:163)  use  the  simplifying  assumptions  that  T[y)  =  T{\)  and  Dxs  is 
independent  of  x.  Both  use  the  assumption  that  xl   <2C  1. 

The  present  study  requires  that  x  be  allowed  to  vary  from  0  to  xe  (up  to  near  1)  in 
order  to  adequately  examine  the  AFWL  trapped  electron  data  base  (Ref.  37).  In  addition, 
if  x  is  not  small,  the  approximation  that  T(y)  =  T{1)  is  a  poor  one,  since  y  is  not  close  to 
1. 

The  function  T{y)  can  be  shown,  within  0.57??  to  be  approximated  as  (Rrf.  10:4030) 

(11)  T[y)  a,  r(0)-[T(0)-r(l)]y'V4 

where  T(0)  =  1  +  Uln(2  +  v/3))/(2n^)    W  1.3802  and  T{\)  =  (srV^Vs  Pa  0.7405  (Ref. 
45:19). 

Obviously,  if  x  approaches  1,  y  approaches  0  and  the  T{y)  =  T(l)  approximation 
becomes  invalid. 

The  three  limiting  approximations  are  removed  in  the  treatment  of  Shuli  (Ref.  47)  by 
introduction  of  a  new  "canonical"  variable  z  such  that: 

(12)  z  =  Z{y)  =    /   x'T{y')dx'  =    /   y'  T{y')  dy' 

and  a  corresponding  diffusion  coefficient 
,13)  Dzz  -  [xT(y)]-DIx. 
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Figure  2.  Postulated  Decay  cf  Higher  Flux  Eigenmodes  to  the  Steady  Strtc  by  Pitch- 
Angle  Diffusion  (Ref.  41:313). 
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Figure  3.  Temporal  evolution  of  1.0  Mev  F,l<y;tronn  from  the  Russian  3  Burst  fRef. 
41:310). 
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Using  (11)  in  (12)  and  performing  the  integration: 


z  «  J  /[r(o)  -  [r(o)  -  m)]/'* 


(14)  z  « 


~  J  \T{0)y'-[T{<S)-T{\)}/' 

r(o)(i-y2)     4[r(o)-ni)][i-y11/4] 


2  11 

whore  7(0)  =  18/35  and  Z{\)  =  0  (Ref.  46:5213). 

Using  (7),  (12)  and  (13),  the  diffusion  equation  becomes: 


2Z  -  i.fD   °J 
(15)  at  ~  a*L      a* 


+  5. 


Shah  (Ref.  47:8)  states  that  (15)  is  a  "canonical"  diffusion  equation  in  that  there  is  no 
Jacobian  factor  which  fails  to  commute  with  Dzx  as  was  the  case  in  (7). 

The  assumptions  which  now  must  be  made  are: 

(1)  Dzm  is  some  "suitably  simple  function  of  z"  (Ref.  47:0)  so  that  exact  ei  nonfunctions 
gn(z)  may  be  specified, 

(2)  5,  the  distributed  natural  environment  source,  is  independent  of  z. 

Assumption  (1)  is  the  limiting  assumption,  since  if  the  function  of  z  is  very  :omplex,  the 
diffusion  equation  becomes  extremely  difficult  to  solve.  Assumption  (2)  is  both  simple  and 
reasonable,  however.  The  source  of  the  natural  environment  must  be  close  tc  constant  in 
order  to  be  the  "driver"  of  quiet-time  equilibrium.  It  must  also  be  distributed  'airly  closely 
to  the  steady-state  quiet  time  distribution  for  the  same  reason. 

Shulz  (Ref.  47:7)  states  that  even  if  Dxx  is  not  of  a  functional  form  tc  yield  exact 
eigenfunctions,  there  may  be  a  Dzx  which  resembles  Dzz  "closely  enough"  and  for  which 
exact  eigenfunctions  gn(z)  (resembling  gn{z))  are  known. 

The  following  derivation  of  the  exact  eigenfunctions  and  the  corresponding  omnidirec- 
tional differential-energy  fJux  arises  directly  from  Shulz  (Ref.  47). 


The  first  assumption  by  Shulz  is  that 
(16)  I>"asG02>"' 


where  a  is  some  number,  not  necessarily  an  integer,  less  than  2,  and  where  Dxz    is  the  value 
of  Dzx  at  some  ze   <    10/35  where  /  vanishes.  Roberts  (1069)  also  uses  this  assumption. 

This  form  of  Dzz  allows  a  basis  set  of  orthogonal  eigenfunctions  to  be  s!  own  for  the 
interval  of  interest:  0    <    z    <    ze.  Following  the  notation  above,  wo  sock  seme  Dzx  and 


eigenfunctions  ?n(~)  which  "closely  resemble"  the  true  functions  and  -which  sa  isfy 


(17) 


dz 


D, 


drJn(z) 


+  X„3«(*)  =  n 


for 


with  the  substitution 


where 


»— (i)V 

l-<7 


and 


/?  = 


^n~(? 


then  (17)  becomes  Bessel's  equation.  The  exact  eigenfunctions  of  (17)  are  giver,  by 


(18) 


where  Kvn  are  the  nth  zeroes  of  the  Bessel  function  Ju  of  the  order  i>,  where 

a  -  1 


u  = 


2-<r' 


The  yn  given  by  (18)  are  normalized  so  that 
(18)  a 


f        On  9m  dz  =  6nm. 
JO 


The  eigenvalues  of  (18)  arc  given  by 


Shulz  and  Boucher  (Ref.  48:fi)  state  that,  since  particles  are  not  lost  at  z  =  0  (but 
rather  at  z  =  rc),  the  diffusion  current  must  go  to  zero  at  z  =  0.  A  diffusion  current  at 
2  =  0  would  imply  diffusion  of  pitch  angles  (and  hence  mirror  points)  into  '-he  equator, 
which  is  not  reasonable.  The  condition  of  no  diffusion  current  at  z  =  0  co 'responds  to 
liraz_0  Ds,  g'n(z)  =  0. 

Shulz  and  Boucher  (Ref.  48)  use  the  series  expansion  of  Jv  from  Abramowic  ;  and  Stogun 
(Ref.  1)  to  show  the  limit  as  z  approaches  zero.  This  author  shows  the  limit  directly  for  the 
special  case  eigenfunction  {J—\/*)  which  is  used  in  this  report.  The  derivatior  is  shown  in 
Appendix  A. 
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In  order  to  arrive  at  the  differential-energy  omnidirectional  flux,  Shulz  (Ref  47)  assumes 
that  f(z,  t)— *fao{z)  as  t— »oo;  that  is,  the  solution  distribution  function  has  some  steady-state 
value.  This  also  assumes  that  5  is  constant  in  time,  but  not  in  x  or  z.  The  f^  distribution 
is  necessary  in  order  to  expand  7  as  a  series  of  eigenfunctions  in  the  flux  equation.  To  arrive 
at  an  expression  for  7oo  we  start  by  setting  Of/Ot  to  zero  in  (15)  and  integrating  twice, 
which  gives: 

(2°)  ?»W-  /"tt-T  S(*")  dz"dz>. 

Jz      Uz'z'   JO 

Expanding  5'  as  a  series  of  orthogonal  eigenfunctions  gives: 

(20a)  S(z")=  f>n  ?„(/') 

n— 0 

and  using  orthogonality  (18a)  with  (20a)  gives: 

(20b)  am  =  fo's{z"')Tjm{z"')dz'". 

Substituting  (20a)  into  (20)  gives 

(20C)  iM-PThL      S   «nTJn(z")dz"dz>. 

Jt    DssJo    ^ 
Rearranging  the  order  of  summation  and  integration  in  (20c)  gives: 

(20d)  7oo(*)=f>»/'-^-/   dn(^")dz"dz'. 

Now  integrating  (17a)  twice  gives: 

or  since  £n  is  a  constant, 

(20e)  m.fw-^tJf)^^. 

Substitution  of  (20e)  into  (20d)  gives: 

(20f>  7«w-E«-!!H- 

n-0  Af» 

Substitution  of  (20b)  into  (20f)  gives: 

(  21 )  U{z)  =  f;  l^±rs(z')  !.(/)  dz> 

„-o     X«    •/o 

where  ${*)  is  positive  (hence  a  source). 

Shulz  and  Lanzerotti  (Ref.  45:39-40)  point  out  that  the  directional  flux  j(E;r)  = 
P~  f\P\tV±.\  ?)  where  /  is  the  same  distribution  function  as  in  this  paper.  SI)  ilz  (Ref.  47) 
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points  out  that  j  must  be  integrated  over  all  angles  in  the  unit  sphere  in  p-sprce  to  obtain 
the  omnidirectional  flux  J: 


(22) 
Since 

and 

it  follows  that 

and  that 
Now  since 


/•cos  a, 

J  ■=  4xp2  /  /  (f(cos  a 

Jo 


sin  a       sin'ao 


B 


BD 


cos 

2. 


.  ,   (B0\ 

cosa  =  y'l-r(|) 

-l(i?/Z?oMy2) 
alcoB  a)  =  — ;        —=^±=±- . 
\fT^j-(B/B0) 


l=I-|(l--)  =  1-|^ 


if  cos  a  =  cosae,    y  =  ye2;  if  cos  a  =  0,    y  =  Z?o/5  and  (22)  becomes 

-Hfl/ao)7d(y2) 


'B)    yJX-y^B/Bo) 
which,  upon  combining  terms  and  reversing  the  limits  (-): 

(23)  j  =  ^,[(B"IB){BJ£^2Ml 

y/l-y-(B~/B0)' 


Integrating  (23)  by  parts 


I  udv  =  uv  —  I  v  du     where  u  =  J 
(B/B0)d{y2) 


dv  = 


y/l-yHB/Bo)' 


gives: 


2?rp- 


v  =  -2\   \ 


{B,/B) 


-<!) 


and  since 


d~J  0~J     dz 


d(y*)        dz  ,%-') 
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■will  give  (since  the  definite  integral  evaluates  to  zero,  because  /  is  0  at  ye  by  definition): 
(24) 


»[iB'/B\  [T    JBn\df     dz     .,« 


J  =  +2jrp 

But,  from  (12),  dz  =  -y  T(y)  dy.  So 

dz  -yT{y)dy        -yT{y)dy        -T(y) 


d(y-)  <%- 

Hence,  (24)  ultimately  becomes: 
(25) 


2y  dy 


— '"C'V-K!)™^- 


Shulz  then  expands  f  xs  a  weighted  series  of  eigenfunctions: 
(2ft)  f  =  f<30{z)+YtAn(E,L;t)yn(z). 


n— 0 


This  allows  the  flux  J  in  (25)  to  be  written: 
(27) 


/  =  2tp- 


r{Bo/D)      /—  f   R  . 


-Y,An(E,L,t)J^         Ji-^ILJ  T{y)?n(z)d(y>) 


The  following  two  expressions  are  derived  in  Appendix  B 
(28) 


l-(»/2) 


n-0     X"       ^ 
0)  

*(;) ,4-"(r'"fcc-;)(£) 


The  actual  form  of  the  flux  ecjuation  awaits  the  choice  of  a,  which  also  determines 
the  dependence  of  Dsz  on  z.  Shulz  uses  a  =  0  for  purposes  of  comparison  of  several 
computational  methods.  This  choice  is  also  attractive  because  it  makes  DXM  independent 
of  z  (see  17(b))  which  was  a  major  assumption  to  begin  with: 


(30) 


D;;=D„t. 


This  choice  is  made  in  the  present  study  for  the  above  reasons  and  also  for  the  sake  of 
simplicity  in  computations. 


With  <t  =  0,  the  Bessel  functios  within  eigenfunctions  become  order  —  1/ !,  whkh  can 
be  simply  written  as: 


(31) 

(32) 


J-iMR)  =  y^R  coaW' 


/'_1/2(i2)  =  -\/  — sin(ii) 

V    TTxl 

and  the  zeroes  /cn  occur  at  (2n  +  1)(t/2). 


cos  [R) 


nR-sfll{xR) 


This  simplifies  (32),  since  the  second  term  will  always  be  zero  if  J'_  w2(k„)  is  evaluated 
at  Kn  =  (2n  +  1)(jt/2). 


Equation  (32)  becomes 

-/_I/i(«n)=  - 


?r(2n+  l)(?r/2) 


sin((2n  +  l)|) 


(33) 


/LiW= 


>r(2n  +  l)(x/2) 


(-D 


n+l 


The  integral  with  respect  to  /  within  (28)  would  thus  be  evaluated  as: 


y«.  y/2/^  y/?7^  >/,(flw.,lH.8/SM.TO  A,    An    *       V 

/      . cos    (2n  +  1)  — -  y    rfr' 

=    /"•  v^  cos((2n  +  l)(T/2re)  ^)dr' 
jo 


(34) 


-sin((2n+  1)(t/2)) 
^27^3ia((2n-H)(ff/2)(z/^)) 

_    y/2j7e  sin((2n+l)(x/2)) 

n+l    • 


((2n+l)(T/22e))(- 

The  Xn  in  (28)  would  be  evaluated  as  (from  (19)): 


(35) 


Xn  =  (2n+l)-(^-)  D„. 


The  flux  equation  (27),  using  (31),  (33),  (34),  and  (35),  will  then  simplify  to  the  following 
for  a  =■  0: 
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(30) 


J  =  2tP: 


(B9/U)  Sy/[2/zt)  oxnUln+mw/'i)) 

E«  (2n+l)|ir/2i,)(-l)"^' 

n— 0 


(2n+L)"(Jr/22er^„, 


V^ 


V*(2n+U(*/a)  (    lJ 


t{(b:X*)  'v^TTFTiK^)™5!'2"-'1'-^) 


(~  vU» + i)(U)W^  *(•*"  * 1}  i  *) 


CM((8n+l)(x/2)(z/*„)) 


«««»+lK*/«X*/«.))*  v/.,.^,,,8,/,.!,,/,.! 


JLfiY"  / 


1)(t/2)(3/^ 


COS 


(*♦' >§;) 


+  ^(_sin((2n+1)ii)^Z 


+  l)(r/8K*/*) 


co»((2n+l)Qr/2)(g/*e)) \ 

n  +  l)(T/2)(z/*e))V2/(T(2n+  l){z/ze))J 


(- 


Equation  (36)  reduces  to  the  following  upon  simplification  and  cancellation  of  terms 
(37) 
J 


-Art  ?«ss?  •  \RD  tw-w" 


•\/i 


-»•(£)  n»)<rt 
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The  first  term  in  Equation  (37)  represents  a  steady  state  value  of  the  flux,  and  the  second 
term  represents  a  sum  of  eigenmodes.  Each  A„  (Ref.  45:102)  has  the  form  an(/')e~  n'B"- 
As  n  increases,  the  r„'s  become  shorter;  i.e.,  the  higher  eignmodes  decay  fatter,  so  that 
eventually  only  the  fundamental  mode  is  left,  which  decays  exponentially.  At  t  =  oo,  the 
fundamental  mode  has  reached  zero,  and  only  the  flux  remaining  is  the  steady  state  flux, 
which  corresponds  to  the  natural  environment  background. 

The  remaining  substitutions  in  (37)  are  ze,  z,  T(y),  and  ye: 


(38)  z  =  Z(V)  =  i(l  -  r)T[0)  -  1[T(0)  -  T(l)](l  -  y<ll/<>), 

zc  =  *(*)  -  5(1  -  !/,2)r(o)  -  l[7i(o)  -  r(i)](i  -  *<»/«)) 

(as  shown  in  Equation  (14))  T(y)  is  given  by  Equation  (11). 

The  variable  of  integration  in  (37)  is  taken  as  y  rather  than  y2  for  simplicity  of 
integration. 

The  preceeding  derivation  has  explicitly  followed  Shulz  (1981).  However,  Shulz  did 
not  show  how  the  differential-energy  flux  related  to  the  integral-energy  flu>  .  Shulz  and 
Lanzerotti  (Ref  45:163)  state  that  the  integral  flux  "will  then  scale  as"  the  diff?rential  flux, 
without  derivation.  Since  the  AFWL  trapped  electron  data  base  (Ref.  37)  exists  primarily  as 
omnidirectional  Duxes  integrated  over  broad  energy  bandpass,  the  relation  must  bo  derived. 

It  will  be  shown  in  the  following  and  Lanzcrotti  (Ref.  45)  are  indeed  correct:  The 
differential-energy  flux  equation  (37)  will  scale  as  the  integral  flux,  and  the  energj  dependence 
will  be  imbedded  in  the  constants  for  the  source  and  the  eigenmodes. 

From  Roederer  (Ref.  39:86),  the  integral  energy  omnidirectional  Dux  can  be  written: 
(39) 


POO 


J>E=  I    JdE 

Je 


where  /  is  the  differential-energy  flux,  as  in  the  previous  equations,  and  J>e  is  the  integral 
energy  flux  above  threshold  energy  E. 

A  principal  assumption  will  be  that  energy  and  pitch  angle  are  both  independent 
variables  of  the  flux.  Thus,  (26)  would  really  look  like 

(39  a)  j{ZtE)=z-fooi{E)  J^  +   jrAn{E,L;t)TJn(z) 


n  —  0 


and  hence  in  (25)  Of/dz  would  look  like 

(39b)  dl  =  -Iooi{E)  7L.W-+    £  An(E,L;t)g'n(=). 

n—Q 


Using  (27),  (39a)  and  (39b),  Equation  (39)  becomes,  at  a  given  //-shell: 


(40) 

Je  Jz,  ^0Je  J,t 

where  7oo,(x)  aD(^  5n(z)  represent  ill  of  the  angular-dependent  functions  of  x,  v,  or  z  shown 
in  (25)  through  (27). 

Using  the  previously  stated  assumption  (Ref.  45:1(52-103)  that  An, \E)  are  of  the  form 
an{E)e-t^T"[h:)) ,  40)  becomes: 

J>E  =  -2*fE    i-~Ux{E)dE  f^l^(z)dz 
(41) 

—    it  e 

- 

at  a  given  L-shcll. 


-  2*f;   r  jfianWe-t^W  dE  /'%'„(*)<** 

„  J  E  J 1 1 


To  "remove"  the  energy  dependence  of  the  decay-time  terms,  the  expectation  value 
of  e~ ,/'r")  is  computed,  assuming  that  the  second  term  in  (-11)  represents  a  listribution- 
function  of  the  energies. 

(42)  V  /  fiVM*)« 

Then  (41)  can  be  written,  using  (42): 
J>E  =  1x\[-jE    ?7mi{E)4B[    jz  *~Lt(x)dz 

(«)  -2*T,[f~ P*«n{E)dE      L-*l^     f\n{x)dT. 

Now  since  it  is  assumed  that  the  energy  and  angular  shape  functions  are  separable, 
the  foo^iE)  must  include  only  those  parts  of  (28)  which  have  energy  dependence.  The  only 
terms  in  (28)  which  could  have  any  energy  dependence  are  5  and  X„. 

Shulz  (Ref.  47:23)  notes  that,  while  /^(z)  resembles  jo(-)  in  functional  form,  / «(~) 
will  coincide  exactly  with  go(z)  only  if  3  is  directly  proportional  to  30(2)-  Per  simplicity, 
the  assumption  is  made  in  this  study  that  5  is  constant  over  the  interval  0  <  2  <  zc.  In 
physical  terms,  this  really  corresponds  to  a  steady-state  isotropic  source.  It  is  known  that 
the  natural  environment  is  not  isotropic  over  either  space  or  time,  but  the  viriations  are 
generally  small  compared  to  the  variations  from  an  injection  by  a  nuclear  evei  t. 

Now  if  5  and  X„  are  the  only  terms  in  (28)  which  may  have  energy  depe  idence,  then 
foe>l{E)  in  (43)  would  correspond  to 

f*  SjS) 

and  since  the  only  portion  of  the  Xn  which  may  have  any  energy  dependence  ("rom  (10))  is 
DI:,,  then  7eoi(^)  should  be  proportional  to 

S(E) 
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where  DtSt(E)  may  be  a  constant  in  E.  Then?  0137,  of  course,  be  some  other  functional 
energy-dependence  included  in  focl(E)  which  is  yet-undetermined. 

Equation  (43)  is  generally  of  the  same  functional  form  as  equation  (37)  if  Lhc  terms  in 
square  brackets  in  (43)  correspond  to  the  combined  constant  terms  in  (37)  as  demonstrated 
above  and  as  shown  in  Appendix  D.  Equation  (13)  also  bears  out  Shulz  and  Lanzjrotti's  (fief. 
45:183)  prediction  that  the  integral  flux  scales  with  the  differential  flux;  i.e.,  the  "pitch-angle 
shape"  dependence  is  invariant  for  differential-  or  integral-energy  fluxes. 

Since  it  is  impractical  to  evaluate  an  infinite  number  of  terms  in  (37),  only  the  terms 
up  to  n  =  3  will  be  shown.  Shulz  and  Lanzerotti  (Ref  45:1(33)  state  that  the  I  igher  modes 
vnnish  for  n  >  2.  This  in  fact  means  that  the  higher  modes  decay  so  rapidly  that  modes 
higher  that  n  =  3  should  not  be  seen  except  at  very  early  times.  The  final  form  of  liquation 
(43),  taken  to  n  =  3,  is  shown  in  Appendix  C.  This  is  the  form  of  the  flux  equ:.tion  used  in 
the  data  analysis,  except  that  fewer  terms  were  used  than  are  shown  in  Apperdix  C. 

The  general  solution  for  the  differential-energy  flux  (which,  from  (43),  is  directly  propor- 
tional to  the  integral-energy  flux),  using  Bessel  functions  of  arbitrary  order  and  with  a- 
dependence  explicitly  shown,  is  given  in  Appendix  D. 


IV.  ANALYSIS  OF  DATA 

The  AFWL  Trapped  Electron  Data  Base  consists  of  a  set  of  computer  tapes  of  satellite 
counting  data  which  have  been  assembled  from  historical  archives  (such  as  I  he  National 
Space  Sciences  Data  Center,  NSSDC)  (Ref.  37).  The  data  cover  injections  of  electrons  from 
nine  high- altitude  nuclear  detonations  shown  in  Table  1  (Ref.  6:6-2). 

The  Data  Base  has  been  organized  by  Pfitzer  into  a  coherent  set  of  tapes  which 
present,  by  satellite,  net  electron  omnidirectional  integral-energy  (above  threshold)  fluxes  as 
a  function  of  time,  B,  and  L-shell.  The  computed  errors  for  each  data  point,  in  units  of  flux, 
are  also  given.  The  background  which  was  subtracted  to  give  net  flux  is  also  given  (Ref. 
37:11,  188-189). 

Additional  data  have  been  collected  by  AFWL  from  plots  published  in  early  literature. 
These  plots  have  been  photographically  enlarged,  and  the  data  points  digitized  and  placed 
in  computer  files  with  formats  similar  to  the  Pfitzer  data  (Ref.  12:670-671,  Ref.  31:646-468, 
Ref.  38:637-638). 

The  satellites  which  provided  the  data,  and  their  orbital  parameters,  are  shown  in  Table 
II  (Ref.  44:41-45). 

The  analysis  of  the  data  base  required  development  of  three  computer  programs  by 
personnel  at  AFWL,  in  coordination  with  this  author: 

(1)  "Program  DTABASE"  to  read  raw  data,  organize  it  for  processing,  aad  output  it 
into  uniformly  formatted  files  (Ref.  15). 

(2)  A  program  to  plot  the  raw  data  and  to  plot  flux  curves  generated  by  the  fitting 

program. 

(3)  "Program  Electrofit",  a  fitting  program  to  take  data  points  from  DTVBASE  at  a 
given  L-shell  and  energy  3nd  to  fit  them  with  the  theoretical  model  leveloped  in 
this  study.  It  also  generates  flux  curves  using  the  fitted  functions  at  predetermined 
times  to  compare  with  the  raw  data  (Ref.  34). 

The  first  two  programs  were  written  with  limited  input  from  the  author.  T'le  third  wns 
written  under  close  and  extensive  coordination  with  this  author,  as  discussed  in  this  chapter. 
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The  DTABASE  program,  written  by  Miss  Cherise  Jarrett  of  Computer  Sciences 
Corporation,  reads  the  raw  data  Qles,  in  various  formats  created  by  PBtter  (Rcf.  37),  satel- 
lite by  satellite,  and  stores  them  in  a  three-dimensional  array.  The  points  are  Lhen  ordered 
(within  each  energy  group)  by  increasing  time  since  burst.  To  generate  a  file  of  points  at  a 
specified  L-value,  each  data  point  lvalue  is  compared  to  the  desired  L  and  is  recorded  in 
the  output  file  if  a  match  occurs.  Additional  points  are  "created"  by  linear  interpolation  of 
two  successive  data  points  if  they  fall  on  both  sides  of  the  desired  L  and  if  the  data  points 
differ  in  time  by  no  more  than  a  specified  value.  The  interpolation  is  perforned  in  L,  B, 
time,  flux,  and  error.  The  output  file  can  be  restricted  to  a  speciGed  time  "window"  and  to 
a  specified  range  of  B  values,  if  desired;  however,  the  data  were  analysed  over  the  full  range 
of  B  and  over  all  times  which  existed  in  the  Data  Base.  The  program  has  tie  additional 
capability  to  "/^-average"  data  points  at  the  same  L  if  their  individual  B  vahn  s  are  within 
a  user  specified  limit  of  each  other  and  they  are  within  the  required  time  difference.  This 
procedure  has  the  effect  of  eliminating  "double"  points,  and  was  used  for  the  majority  of 
this  study  with  a  limit  of  one  percent.  Another  capability  of  the  program  is  to  interpolate 
to  a  specified  0-value  in  a  manner  analogous  to  the  //-interpolation  described  above;  this 
capability  was  unused  in  this  study  (Ref.  15:1-2). 
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TABLE  n 


Satellites  and  Their  Orbital  Parameters 


Name  Launch      Failure       Decay       Period     Perigee     Apogee      Incl. 

Date  Date  Date         (rain.)        (mi.)  (mi)       (^e?) 


Alouette  1 

9-28-62 

unk. 

— 

105.4 

020 

63? 

80.5 

Alpha 

Upsilon  1 

9-1-62 

— 

10-26-64 

94.4 

189 

41? 

82.8 

Explorer  15 

10-27-02 

2-9-63 

— 

312.0 

194 

10710 

18.0 

Injun  1 

6-29-61 

3-6-83 

— 

103.8 

534 

631 

67.0 

Traac 

11-15-61 

7-62 

— 

105.6 

562 

720 

32.4 

Telstar 

7-10-62 

2-21-63 

— 

157.8 

593 

35C3 

44.8 

Explorer  4 

7-26-58 

10-8-58 

10-23-59 

110.1 

163 

(Ref. 

1372 
44:41-45) 

50.1 
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In  order  to  determine  the  optimum  time  difference  to  use  for  interpolation  across  the 
desired  L- value,  the  DTABASE  program  was  modified  under  direction  of  this  author.  The 
program  now  presents  a  summary  table  of  the  raw  data  points  for  each  satellite,  showing  the 
frequency  of  occurance  of  time  differences  between  successive  data  points  from  0.1  minute 
to  5.0  minutes,  by  tenths  of  a  minutes.  The  cumulative  percentage  of  points  with  a  specified 
time  difference  or  less  is  also  shown.  These  summaries  are  presented  for  the  bursts  and 
satellites  examined  as  Tables  111  through  XIII  (Ref.  15:67-70).  Perusal  of  those  t.ibles  showed 
that  the  optimum  time  difference  for  interpolation  was  2.2  minutes,  primarily  to  include  the 
maximum  number  of  points  in  the  Starfish  Telstar  file.  This  time  difference  was  used  for  all 
data  examined  in  this  study. 

A  final  capability  of  DTABASE  is  to  add  data  points  to  a  specified  Rvalue  file  by 
rounding  points  in  L  to  the  desired  value.  This  procedure  assumes  that  poin.s  which  are 
closer  to  the  desired  L  than  some  specified  limit  may  be  considered  to  occur  a!  the  desired 
L.  Within  the  errors  in  satellite  measurements,  this  procedure  is  reasonable,  and  in  fact 
appears  (by  data  comparison)  to  have  been  used  by  Roberts  (Ref.  41:310),  although  he  does 
not  so  explicitly  state.  This  capability  exists  in  the  program,  but  was  not  used  n  this  study 
(Ref.  15). 

The  plotting  program  was  developed  by  Mr.  John  Burgio  of  AFWL  to  generate  all  the 
plots  of  raw  data  and  fitted  curves  shown  in  this  report.  The  program  plots  raw  data  points 
beginning  at  a  specified  start  time  since  burst  (in  days)  and  covering  the  time  period  in  a 
specified  time  "Svindow"  (shown  in  Figure  headings  as  "TW").  The  lvalue,  equatorial  B- 
value,  and  B-cutoff  values  are  also  shown.  If  the  /?-range  of  points  is  restricted,  the  minimum 
and  maximim  values  of  B  are  shown.  The  equatorial  B  value  is  computed  from  L  by  the 

relati°a  „         0.308 

*--zr- 

The  tf-cutoff  value  was  initially  computed  by  using  the  dipole  formula  (Ref.  3?):55) 


1 


LV4-(3//,) 


-  y/{Bo/Beut). 


However,  it  was  later  determined  that  this  approximation  might  be  less  than  optimum;  so 
a  table  of  L-values  versus  Beut  was  used  and  was  linearly  interpolated  across  the  tabulated 
//-values  to  find  the  Beut  value.  This  table  (Table  XIV)  was  prepared  by  AFWL  from  the 
48-term  Jensen-Cain  model  of  the  magnetosphere  (Ref.  16).  It  was  found  that  use  of  Table 
XIV  produced  Beut  (or  yeut  or  zeut)  values  more  indicative  of  the  global  cutoff  values.  This 
table  was  used  thereafter  in  both  the  plotting  program  and  the  fitting  program  (discussed 
below).  This  procedure  is  reasonable  in  view  of  the  use  of  the  same  Jenscn-C  lin  model  to 
produce  the  B-L  coordinates  for  the  PfiUer  Data  Base  (Ref.  37). 


The  plotting  program  plots  flux  as  a  function  of  z  =  \/l  —  (Bo/B).  It  also  has  the 
capability  to  plot  flux  data  points  versus  time  over  all  values  of  L  and  B,  or  flux  versus 
energy  at  specified  B  or  L.  Next  to  the  satellite  name  on  each  plot  which  has  fitted  curves 
is  a  number  in  wjnar"  bracket?  which  corresponds  to  the  serial  number  of  the  run  of  the 
fitting  program  which  prodeuced  the  curves. 

This  author's  input  to  the  plotting  program  was  limited  to  supplying  header  information 
and  specifying  layout  and  format. 

The  third  program  used  in  data  analysis  is  "Program  Elcctrofit,"  written  by  Mr.  narry 
Murphy  of  AFWL  with  the  close  coordination  of  this  author  (Ref.  34).  Electro.1t  is  a  least- 
squares  fitting  program  which  uses  the  functions  developed  by  this  auLlmr  (Appendix  C) 


to  fit  data  points  which  are  output  from  DTABASE  in  order  to  compute  the  characteristic 
decay  times  (r„)  and  linear  amplitude  constants  (a„).  These  decay  limes  and  constants  are 
determined  at  the  specific  energy  and  Rvalue  which  are  input  with  the  data. 
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Table  HT 


COUNTS    OF    TIME    DIFFERENCES    BETWEEN   NONZERO   RAW    DATA    POINTS 
FOR    ALL    CHANNELS   OF    SATELL ITE < S ) :        TELSTAR 
FOR    BURST         STARFISH  p 

NO 


TIME 

BIN 

(MIN) 

0    0 

0    1 

0    1 

0    2 

0    2 

0    3 

0    3 

0    4 

0    4 

0    5 

0    5 

0    6 

0    6 

0    7 

0    7 

0    8 

0    8 

0    9 

0    9 

1    0 

1    0 

1     1 

1     1 

1    2 

1.  2 

1     3 

1.  3 

1.  4 

1.  4 

1.  S 

1.  5 

1    6 

1.  6 

1.  7 

1    7 

1.  8 

1.  S 

1    9 

1.  9 

2    0 

2    0 

2.  1 

2    1 

2.  2 

2    2 

2    3 

2.  3 

2    4 

2    4 

2    5 

2    5 

2.  6 

2.  6 

2.  7 

2.  7 

2    8 

2.  8 

2    9 

2    9 

3.  0 

3    0 

3    1 

3.  1 

3.  2 

3    2 

3    3 

3    3 

3    4 

3    4 

3    5 

3    5 

3.  6 

3    6 

3    7 

3.  7 

3.  8 

3    8 

3    9 

3    9 

4.  0 

4    0 

4    1 

4    1 

4.  2 

4    2 

4    3 

4    3 

4    4 

4    4 

4    5 

4    5 

4    6 

4    6 

4    7 

4.  7 

4    8 

4    8 

4    9 

4    9 

5.  0 

5   OO 


DELTA    TIMES 

CUMULATIVE 

IN    BIN 

PERCENTAGE 

607 

1     4 

3 

1     4 

7 

1     5 

33 

1.  5 

0 

1.  5 

22 

1     6 

9 

1    6 

58 

1.  7 

2 

1     7 

219 

2    3 

29 

2    3 

4 

2    3 

56 

2.  5 

12 

2    5 

15 

2    5 

0 

2    5 

35 

2    6 

11 

2    6 

19 

2    7 

33421 

81.  4 

178 

81     8 

26 

81     9 

8 

81     9 

5 

81     9 

1 

81     9 

5 

81     9 

0 

81.  9 

2 

81     9 

2 

81     9 

50 

82    0 

6 

82     1 

0 

82    1 

3 

82    1 

1 

82     1 

1 

82     1 

0 

82     1 

2 

82     1 

1 

82.   1 

5 

82    1 

1451 

85    5 

5 

85    5 

3 

85    5 

0 

85    5 

0 

85    5 

0 

85    5 

2 

85    5 

6 

85    5 

'2 

85    5 

2 

85    5 

23 

85    6 

61  15 

100    0 

20 


T.ible  IV 


COUNTS   OF    TIME   DIFFERENCES   BETWEEN   NONZERO   RAW   DATA    POINTS 
FOR    ALL    CHANNELS   OF    SATELL ITE( S) :        EXPLORER4 
FOR    BURST         ARCUS    1.2     ,3 


TIME 

BIN 

NO       DELTA   TIMES 

CUMULATIVE 

(MIN) 

IN    BIN 

PERCENT AOE 

0    0 

0    1 

442 

63    7 

0     1 

0    2 

133 

82    9 

0    2 

0    3 

46 

89    5 

0    3 

0    4 

10 

90    9 

0    4 

0    5 

8 

92     1 

0.  5 

0    6 

0 

92    1 

0.  6 

0    7 

1 

92    2 

0.  7 

0    8 

1 

92    4 

0    8 

0    9 

0 

92    4 

0    9 

1    0 

0 

92    4 

1     0 

1     1 

0 

92    4 

1     1 

1    2 

0 

92    4 

1    2 

1.  3 

0 

92    4 

1    3 

1     4 

0 

92    4 

1    4 

1     5 

0 

92    4 

1     5 

1.  6 

0 

92    4 

1    6 

1.  7 

0 

92    4 

1.  7 

1    8 

0 

92    4 

1.  S 

1.  9 

0 

92    4 

1.  9 

2    0 

0 

92    4 

2    0 

2    1 

0 

92    4 

2    1 

2.  2 

0 

92.  4 

2    2 

2    3 

0 

92    4 

2    3 

2    4 

1 

92    5 

2.  4 

2    5 

0 

92    5 

2    5 

2    6 

0 

92    5 

2    6 

2    7 

0 

92    5 

2    7 

2.  8 

1 

92    7 

2    S 

2    9 

1 

92    8 

2    9 

3    0 

0 

92    8 

3    0 

3    1 

1 

92    9 

3    1 

3    2 

1 

93    1 

3    2 

3    3 

0 

93    1 

3    3 

3    4 

0 

93    1 

3    4 

3    5 

0 

93     1 

3    5 

3    6 

0 

93    1 

3    6 

3    7 

0 

93     1 

3    7 

3    8 

0 

93    1 

3    8 

3    9 

0 

93     1 

3.  9 

4    0 

0 

93    1 

4    0 

4.  1 

0 

93    1 

4     1 

4    2 

0 

93    1 

4    2 

4    3 

0 

93    1 

4    3 

4    4 

0 

93    1 

4    4 

4    5 

0 

93     1 

4    5 

4    6 

0 

93    1 

4    6 

4    7 

0 

93    1 

4    7 

4    8 

0 

93    1 

4    8 

4    9 

0 

93    1 

4.  9 

5    0 

0 

93    1 

5. 

0<= 

48 

100    0 

~f> 


Table  V 


COUNTS    OF    TIME    DIFFERENCES    BETWEEN    NONZERO   RAW    DATA    POINTS 
FOR    ALL    CHANNELS    OF    SATELL I TE ( S ) :        ALLOUETTE1 
FOR    BURST  RUSSIAN     1 

NO. 


TIME 

BIN 

1 

(MIN) 

0 

0 

0. 

1 

0 

1 

0 

2 

0 

2 

0 

3 

0 

3 

0 

4 

0 

4 

0 

5 

0 

5 

0 

6 

0 

6 

0 

7 

0 

7 

0 

B 

0. 

8 

0 

9 

0 

9 

1 

0 

1 

0 

1 

1 

1 

1 

1 

2 

1 

2 

1. 

3 

1. 

3 

1 

4 

1 

4 

1. 

5 

1. 

5 

1. 

6 

1. 

6 

1 

7 

1 

7 

1. 

e 

1 

8 

1 

9 

1 

9 

2. 

0 

2 

0 

2. 

1 

2 

1 

2. 

2 

2 

2 

2. 

3 

2 

3 

2 

4 

2 

4 

2 

5 

2 

3 

2. 

6 

2 

6 

2. 

7 

2 

7 

2. 

8 

2 

8 

2. 

9 

2 

9 

3. 

0 

3 

0 

3 

1 

3 

1 

3 

2 

3 

2 

3 

3 

3 

3 

3 

4 

3 

4 

3 

5 

3 

5 

3 

6 

3 

6 

3 

7 

3 

7 

3 

8 

3 

8 

3 

9 

3 

9 

4 

0 

4 

0 

4 

1 

4 

1 

4 

2 

4 

2 

4 

3 

4 

.  3 

4. 

4 

4 

.  4 

4 

5 

4 

5 

4 

6 

4 

6 

4 

7 

4 

.  7 

4 

8 

4 

8 

4 

9 

4 

9 

5 

0 

5. 

0<« 

DELTA    TIMES 

CUMULATIVE 

IN    BIN 

PERCENTAGE 

0 

0    0 

3413 

72    7 

0 

72    7 

572 

84    8 

180 

88    7 

15 

89    0 

105 

91    2 

0 

91    2 

46 

92    2 

46 

93    2 

1 

93    2 

5 

93    3 

0 

93    3 

15 

93    6 

4 

93    7 

3 

93    8 

5 

93    9 

0 

93    9 

12 

94     1 

1 

94    2 

1 

94    2 

2 

94    2 

0 

94    2 

2 

94    3 

1 

94    3 

0 

94    3 

3 

94    4 

0 

94    4 

0 

94    4 

0 

94.  4 

1 

94    4 

0 

94    4 

0 

94    4 

3 

94    4 

0 

94    4 

0 

94    4 

0 

94    4 

0 

94    4 

1 

94     5 

1 

94     5 

0 

94     5 

3 

94     5 

0 

94     5 

1 

94    6 

0 

94    6 

0 

94    6 

0 

94    6 

0 

94    6 

0 

94    6 

0 

94    6 

255 

100    0 

31 


Table  VI 


COUNTS   OF    TIME    DIFFERENCES    BETWEEN    NONZERO   RAU    DATA    POINTS 
FOR    ALL    CHANNELS   OF    SATELL I TE ( S ) :        EXPLORER15 
FOR    BURST         RUSSIAN    1 

NO. 


TIME 

BIN 

1 

(MIN) 

0 

0 

0 

1 

0 

1 

0 

2 

0 

2 

0 

3 

0 

3 

0 

4 

0 

4 

0 

5 

0 

5 

0 

6 

0 

6 

0 

7 

0 

7 

0 

8 

0 

B 

0 

9 

0 

9 

1 

0 

1 

0 

1 

1 

1 

1 

1 

2 

1 

2 

1 

3 

1 

3 

1 

4 

1 

4 

I. 

5 

1- 

5 

1 

6 

1 

6 

1 

7 

1 

7 

1 

8 

1 

e 

1 

9 

1 

9 

2 

0 

2 

0 

2 

1 

2 

1 

2 

2 

2 

2 

2. 

3 

2 

3 

2 

4 

2 

4 

2 

5 

2. 

5 

2. 

6 

2 

6 

2 

7 

2 

7 

2 

8 

2 

e 

2 

9 

2 

9 

3 

0 

3 

0 

3 

1 

3 

1 

3 

2 

3 

2 

3 

3 

3 

3 

3 

4 

3 

4 

3 

5 

3 

5 

3 

6 

3 

6 

3 

7 

3 

7 

3 

8 

3 

8 

3 

9 

3 

9 

4 

0 

4 

0 

4 

1 

4 

1 

4 

2 

4 

2 

4 

3 

4 

3 

4 

4 

4 

4 

4 

5 

4 

5 

4 

6 

4 

6 

4 

7 

4 

.  7 

4 

8 

4 

8 

4 

9 

4 

9 

5 

0 

5. 

0<= 

DELTA    TIMES 

CUMULATIVE 

IN    BIN 

PERCENTAGE 

0 

0    0 

0 

0    0 

0 

0    0 

2 

1     7 

12 

12    0 

0 

12    0 

5 

16    2 

6 

21     4 

1 

22    2 

0 

22    2 

13 

33    3 

8 

40    2 

0 

40    2 

1 

41     0 

0 

41     0 

0 

41     0 

11 

50    4 

3 

53    0 

0 

53    0 

0 

53    0 

0 

53    0 

3 

55    6 

4 

59    0 

0 

59    0 

0 

59    0 

0 

59    0 

0 

59    0 

0 

59    0 

1 

59    8 

0 

59    8 

2 

61     5 

0 

61     5 

0 

61     5 

3 

64     1 

7 

70    1 

0 

70     1 

3 

72    6 

0 

72    6 

0 

72    6 

1 

73    5 

0 

73    5 

0 

73    5 

0 

73    5 

0 

73    5 

0 

73    5 

1 

74    4 

0 

74    4 

0 

74    4 

0 

74     4 

0 

74     4 

30 

100    0 

Table  VT! 


COUNTS    OF    TIME    DIFFERENCES    BETWEEN    NONZERO   RAW    DATA    POINTS 
FOR    ALL    CHANNELS    OF    SATELL ITE < S ) :        TELSTAR 
FOR    BURST  RUSSIAN     1 

NO. 


TIME 

BIN 

1 

(MIN) 

0 

0 

0 

1 

0 

1 

0 

2 

0 

2 

0 

3 

0 

3 

0 

4 

0 

4 

0 

5 

0 

5 

0 

6 

0 

6 

0 

7 

0 

7 

0 

8 

0 

B 

0 

9 

0 

9 

1 

0 

1. 

0 

1 

1 

1 

1 

1. 

2 

1 

2 

1. 

3 

1 

3 

1. 

4 

1. 

4 

1. 

5 

1 

5 

1. 

6 

1 

6 

1 

7 

1 

7 

1. 

B 

I. 

8 

1 

9 

1 

9 

2. 

0 

2. 

0 

2 

1 

2 

1 

2 

2 

2 

2 

2. 

3 

2 

3 

2 

4 

2 

4 

2. 

5 

2 

5 

2 

6 

2 

6 

2 

7 

2 

7 

2. 

B 

2. 

e 

2. 

9 

2 

9 

3. 

0 

3 

0 

3 

1 

3 

i 

3 

2 

3 

2 

3 

3 

3 

3 

3 

4 

3 

4 

3 

5 

3 

5 

3 

6 

3 

6 

3 

7 

3 

7 

3 

8 

3 

B 

3 

9 

3 

9 

4 

0 

4 

0 

4 

1 

4 

1 

4 

2 

4 

,  2 

4 

3 

4 

3 

4 

4 

4 

4 

4 

5 

4 

5 

4 

6 

4 

6 

4 

7 

4 

.  7 

4 

8 

4 

B 

4 

9 

4 

.  9 

5 

0 

5. 

0<= 

DELTA    TIMES 

CUMULATIVE 

IN    BIN 

PERCENTAGE 

0 

0    0 

0 

0    0 

0 

0    0 

0 

0.  0 

1 

0    0 

0 

0    0 

0 

0    0 

0 

0    0 

0 

0    0 

7 

0.  2 

1 

0    2 

0 

0    2 

0 

0    2 

0 

0    2 

0 

0    2 

0 

0    2 

0 

0    2 

1 

0    2 

2 

0    3 

3637 

86    6 

3 

86    6 

3 

86    7 

0 

86    7 

0 

86    7 

0 

86    7 

0 

86    7 

0 

86    7 

0 

86    7 

b 

86    7 

6 

86    B 

0 

B6    8 

0 

B6    8 

0 

86    B 

0 

86    8 

0 

86    B 

0 

86    8 

0 

86    8 

0 

86    8 

1 

66    9 

134 

90    0 

0 

90    0 

0 

90    0 

0 

90    0 

2 

90    1 

0 

90    1 

0 

90     1 

0 

90     1 

0 

90.   1 

0 

90    1 

5 

90    2 

413 

100    0 

Table  VTJI 


COUNTS   OF    TIME    DIFFERENCES    BETWEEN    NONZERO   RAW    DATA    POINTS 
FOR    ALL    CHANNELS    OF    SATELL I TE ( S ) :         ALLOUETTE1 
FOR    BURST         RUSSIAN    2 

NO 


TIME 

BIK 

1 

<MIN) 

0 

0 

0 

1 

0 

1 

0 

2 

0 

2 

0 

3 

0 

3 

0 

4 

0 

4 

0 

5 

0 

5 

0 

6 

0 

6 

0 

7 

0 

7 

0 

B 

0 

8 

0 

9 

0 

9 

1 

0 

1 

0 

1 

1 

1. 

1 

1 

2 

1 

2 

1 

3 

1 

3 

1. 

4 

1 

4 

1 

5 

1 

5 

1 

6 

1 

6 

1 

7 

1. 

7 

1 

8 

1 

8 

1 

9 

1 

9 

2 

0 

2 

0 

2. 

1 

2 

1 

2 

2 

2 

2 

2. 

3 

2 

3 

2 

4 

2 

4 

2. 

5 

2 

5 

2 

6 

2 

6 

2. 

7 

2 

7 

2 

8 

2 

8 

2 

9 

2 

9 

3 

O 

3 

0 

3 

1 

3 

1 

3 

2 

3 

2 

3 

3 

3 

3 

3 

4 

3 

4 

3 

5 

3 

5 

3 

6 

3 

6 

3 

7 

3 

7 

3 

8 

3 

8 

3 

9 

3 

9 

4 

0 

4 

0 

4 

1 

4 

1 

4 

2 

4 

2 

4 

3 

4 

3 

4 

4 

4 

4 

4 

5 

4 

5 

4 

6 

4 

6 

4 

7 

4 

7 

4 

8 

4 

8 

4 

9 

4 

9 

5 

0 

5 

OO 

DELTA  TIMES 

CUMULATIVE 

IN  BIN 

PERCENTAGE 

45 

1  1 

3107 

78  9 

0 

78  9 

409 

89  1 

91 

91  4 

7 

91  5 

52 

92  8 

0 

92  8 

29 

93  6 

25 

94  2 

3 

94  3 

13 

94  6 

0 

94  6 

7 

94  8 

1 

94  8 

0 

94  8 

6 

94  9 

0 

94  9 

4 

95  0 

1 

95  1 

0 

95.  1 

4 

95  2 

0 

95  2 

1 

95  2 

3 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95.  3 

0 

95  3 

0 

95  3 

0 

95  3 

0 

95  3 

189 

100  0 

Table  IX 


COUNTS  OF  TIME  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELL ITE ( S )  :    EXPLORER15 
FOR  BURST    RUSSIAN  2 


TIME 

BIN 

NO       DELTA    TIMES 

CUMULATIVE 

(MIN) 

IN    BIN 

PERCENTAGE 

0    0 

0     1 

0 

0    0 

0    1 

0    2 

0 

0    0 

0    2 

0    3 

0 

0    0 

0    3 

0    4 

71 

2    9 

0    4 

0    5 

264 

13    8 

0.  5 

0    6 

344 

28    0 

0    6 

0    7 

142 

33    9 

0    7 

0    8 

77 

37    0 

0    8 

0    9 

31 

38    3 

0    9 

1.  0 

83 

41     7 

1.  0 

1     1 

65 

44    4 

1.  1 

1.  2 

68 

47    2 

1    2 

1.  3 

64 

49    9 

1     3 

1.  4 

63 

52    5 

I.  4 

1.  5 

115 

57    2 

1     5 

1.  6 

75 

60    3 

1     6 

1.  7 

125 

65    4 

1     7 

1    8 

36 

66    9 

i   e 

1.  9 

11 

67    4 

1.  9 

2    0 

33 

68    7 

2    0 

2     1 

12 

69    2 

2     1 

2    2 

29 

70    4 

2    2 

2.  3 

32 

71     8 

2    3 

2    4 

17 

72    5 

2    4 

2    5 

42 

74    2 

2.  5 

2    6 

14 

74    8 

2    6 

2    7 

26 

75    8 

2.7 

2    8 

11 

76    3 

2  e 

2    9 

16 

76.  9 

2    9 

3    0 

10 

77.  4 

3    0 

3.  1 

10 

77    8 

3.  1 

3    2 

10 

78    2 

3    2 

3    3 

7 

78    5 

3    3 

3    4 

8 

78    8 

3.  4 

3    5 

22 

79    7 

3    5 

3    6 

23 

80.  7 

3    6 

3    7 

45 

82    5 

3    7 

3.  8 

14 

83    1 

3    8 

3    9 

18 

83    8 

3    9 

4.  0 

27 

84    9 

4    0 

4    1 

15 

85    6 

4     1 

4    2 

10 

86    0 

4    2 

4.  3 

6 

86    2 

4    3 

4    4 

10 

86    6 

4    4 

4    3 

4 

86    8 

4    5 

4    6 

15 

87    4 

4    6 

4    7 

10 

87    8 

4    7 

4    8 

4 

88    0 

4    8 

4    9 

14 

88    6 

4    9 

5    0 

8 

88    9 

5. 

0<= 

269 

100.  0 

35 


Tabic  X 


COUNTS    OF    TIME    DIFFERENCES    BETWEEN   NONZERO   RAW    DATA    POINTS 
FOR    ALL    CHANNELS   OF    SATELL ITE ( S) :        TELSTAR 
FOR    BURST  RUSSIAN    2 


TIME 

BIN 

NO.     DELTA   TIMES 

CUMULATIVE 

(MIN) 

IN    BIN 

PERCENTAGE 

0    0 

0    1 

0 

0    0 

0.  1 

0    2 

0 

0    0 

0.  2 

0    3 

0 

0    0 

0    3 

0    4 

0 

0    0 

0    4 

0    5 

0 

0    0 

0.  5 

0    6 

0 

0    0 

0    6 

0    7 

0 

0    0 

0    7 

0    8 

0 

0    0 

0    8 

0    9 

0 

0    0 

0    9 

1.  0 

2 

0    1 

1     0 

1     1 

0 

0    1 

1.  1 

1.  2 

0 

0    1 

1.  2 

1     3 

0 

0.  1 

1.  3 

1.  4 

0 

0    1 

1     4 

1     5 

0 

0    1 

1.  5 

1.  6 

0 

0    1 

1.  6 

1     7 

0 

0    1 

1     7 

1.  8 

1 

0    1 

i   e 

1.  9 

0 

0    1 

1     9 

2    0 

2244 

85    6 

2    0 

2    1 

49 

87.  5 

2     1 

2    2 

1 

87    5 

2    2 

2    3 

0 

87    5 

2    3 

2    4 

0 

87    5 

2    4 

2    5 

0 

87    5 

2.  5 

2    6 

0 

87    5 

2.  6 

2    7 

0 

87    5 

2.  7 

2    8 

0 

87    S 

2.  B 

2    9 

0 

87.  5 

2    9 

3    0 

3 

87    h 

3.  0 

3.   1 

0 

87    6 

3.  1 

3    2 

0 

87    6 

3    2 

3    3 

0 

87    6 

3    3 

3    4 

0 

87    6 

3    4 

3    5 

0 

87    6 

3.  5 

3    6 

0 

87    6 

3.  6 

3    7 

0 

87    6 

3.  7 

3    8 

0 

87    6 

3.  8 

3    9 

2 

87    7 

3    9 

4    0 

100 

91     5 

4    0 

4.   1 

2 

91     6 

4.  1 

4    2 

0 

91     6 

4    2 

4    3 

0 

91     6 

4    3 

4    4 

0 

91     6 

4    4 

4    5 

0 

91     6 

4    5 

4    6 

0 

91     6 

4    6 

4    7 

0 

91     6 

4    7 

4    8 

0 

91     6 

4    8 

4    9 

0 

91     6 

4    9 

5    0 

0 

91.  6 

5. 

0<= 

221 

100    0 

3fi 


Tnblc  XI 


COUNTS    OF    TIME    DIFFERENCES    BETWEEN   NONZERO    RAW    DATA    POINTS 
FOR    ALL    CHANNELS    OF    SATELL I TE ( S ) :        ALLOUETTE1 
FOR    BURST:        RUSSIAN    3 

NO. 


TIME 

BIN 

(MIN) 

0    0 

0    1 

0    1 

0    2 

0.  2 

0    3 

0    3 

0.  4 

0    4 

0    5 

0    5 

0.  A 

0.  6 

0    7 

0    7 

0    8 

0    8 

0    9 

0.  9 

1    0 

1    0 

1.  1 

1     1 

1.  2 

1.  2 

1    3 

1    3 

1.  4 

1    4 

1.  5 

1    5 

1.  6 

1.  6 

1.  7 

t.  7 

1.  8 

1.  8 

1    9 

1.  9 

2.  0 

2.  0 

2.  1 

2.  1 

2    2 

2.  2 

2    3 

2    3 

2.  4 

2.  4 

2    5 

2    5 

2.  6 

2    6 

2    7 

2    7 

2.  8 

2    8 

2.  9 

2    9 

3.  0 

3    0 

3.  1 

3    1 

3    2 

3    2 

3    3 

3    3 

3    4 

3    4 

3.  5 

3.  5 

3.  6 

3.  6 

3.  7 

3    7 

3    8 

3    8 

3    9 

3    9 

4.  0 

4    0 

4    1 

4.  1 

4    2 

4    2 

4    3 

4.  3 

4    4 

4    4 

4    5 

4    5 

4.  6 

4.  6 

4.  7 

4    7 

4    8 

4    8 

4    9 

4.  9 

5    0 

5. 

0<« 

DELTA   TIMES 

CUMULATIVE 

IN    BIN 

PERCENTAGE 

123 

0    8 

10901 

74    0 

45 

74    3 

1595 

85    0 

491 

88    3 

55 

88    7 

289 

90    6 

11 

90    7 

128 

91     5 

74 

92    0 

8 

92     1 

55 

92    5 

3 

92    5 

62 

92    9 

28 

93    1 

9 

93     1 

18 

93.  3 

0 

93    3 

18 

93    4 

15 

93.  5 

2 

93    5 

10 

93    6 

0 

93    6 

13 

93    7 

7 

93    7 

1 

93    7 

7 

93    8 

0 

93    8 

*4 

93    8 

2 

93    8 

2 

93    8 

5 

93    8 

0 

93    8 

2 

93    9 

0 

93    9 

1 

93    9 

1 

93    9 

0 

93    9 

1 

93.  9 

0 

93    9 

0 

93    9 

0 

93    9 

0 

93    9 

0 

93    9 

0 

93    9 

1 

93    9 

0 

93    9 

0 

93    9 

0 

93    9 

0 

93    9 

911 

100.  0 

37 


Table  Xn 


COUNTS   OF    TIME   DIFFERENCES   BETWEEN   NONZERO   RAW   DATA    POINTS 
FOR    ALL    CHANNELS   OF    SATELL ITE ( S ) :        EXPLORER15 
FOR    BURST:        RUSSIAN    3 

NO. 


Tire 

BIN 

(MIN) 

0    0 

0.  1 

0    1 

0    2 

0.  2 

0    3 

0    3 

0    4 

0    4 

0    5 

0.  5 

0    6 

0    6 

0    7 

0    7 

0    8 

0    8 

0    9 

0.  9 

1.  0 

1.  0 

1     1 

1     1 

1.  2 

1.  2 

1    3 

1.  3 

1    4 

1     4 

1    5 

1    5 

1    6 

1.  6 

1    7 

1    7 

1.  8 

1    8 

1    9 

i.  9 

2    0 

2    0 

2    1 

2.  1 

2    2 

2    2 

2    3 

2.  3 

2    4 

2.  4 

2.  5 

2    5 

2    6 

2    6 

2.  7 

2    7 

2    8 

2    8 

2    9 

2    9 

3    0 

3    0 

3    1 

3    1 

3    2 

3    2 

3    3 

3    3 

3    4 

3    4 

3    5 

3    5 

3    6 

3.  6 

3    7 

3    7 

3    8 

3    8 

3    9 

3    9 

4    0 

4    0 

4    1 

4    1 

4    2 

4    2 

4    3 

4    3 

4    4 

4    4 

4    5 

4    5 

4    6 

4    6 

4    7 

4    7 

4    8 

4    8 

4    9 

4    9 

5.  0 

5. 

0<= 

DELTA    TIMES 

CUMULATIVE 

IN    BIN 

PERCENTAGE 

46 

0    3 

0 

0    3 

77 

0    8 

588 

4    8 

1654 

16    0 

2144 

30    5 

1629 

41     5 

970 

48    0 

274 

49    9 

194 

51     2 

226 

52    7 

210 

54     1 

293 

56     1 

273 

58    0 

474 

61     2 

442 

64     1 

618 

68    3 

313 

70    4 

194 

71     7 

131 

72    6 

82 

73    2 

92 

73    8 

74 

74    3 

54 

74    7 

45 

75    0 

35 

75    2 

53 

75    6 

34 

75    8 

37 

76    0 

40 

76    3 

29 

76    5 

47 

76    8 

52 

77    2 

66 

77    6 

53 

78    0 

79 

78    5 

155 

79    6 

91 

80    2 

127 

81     0 

111 

81     8 

83 

82    3 

117 

83    1 

87 

83    7 

110 

84    5 

73 

85    0 

100 

85    6 

76 

86    2 

63 

86    6 

80 

87     1 

59 

87    5 

1848 

100    0 

.?3 


Table  Xffl 


COUNTS   OF    TIME   DIFFERENCES   BETWEEN   NONZERO   RAW   DATA    POINTS 
FOR    ALL    CHAWELS    OF    SATELL  ITE (  S) :        TELSTAR 
FOR    DURST:        RUSSIAN    3 

NO. 


TIME 

BIN 

(MIN) 

0    0 

0    1 

0    1 

0    2 

0    2 

0    3 

0.  3 

0    4 

0.  4 

0.  5 

0    5 

0    6 

0    6 

0    7 

0.  7 

0    8 

o  e 

0.  9 

0.  9 

1.  0 

1.  0 

1     1 

1     1 

1    2 

1.  2 

1.  3 

1    3 

1.  4 

1.  4 

1.  5 

1.  5 

1.  6 

1.  6 

1.  7 

1.  7 

1.  8 

1.  8 

1.  9 

1.  9 

2.  0 

2    0 

2    1 

2.  1 

2.  2 

2    2 

2.  3 

2    3 

2    4 

2    4 

2.  5 

2.  5 

2    6 

2    6 

2.  7 

2    7 

2    8 

2.  8 

2    9 

2.  9 

3    0 

3    0 

3.  1 

3.  1 

3    2 

3    2 

3.  3 

3    3 

3    4 

3.  4 

3.  5 

3    5 

3.  6 

3    6 

3.  7 

3.  7 

3    8 

3.  8 

3.  9 

3    9 

4    0 

4    0 

4    1 

4    1 

4.  2 

4    2 

4    3 

4    3 

4    4 

4    4 

4    5 

4    5 

4    6 

4    6 

4    7 

4    7 

4.  8 

4    8 

4    9 

4    9 

5.  0 

5 

0<= 

DELTA   TIMES 

CUMULATIVE 

IN    BIN 

PERCENTAGE 

5 

0    0 

0 

0    0 

0 

0    0 

0 

0    0 

0 

0    0 

0 

0    0 

0 

0    0 

0 

0    0 

1 

0    0 

64 

0    3 

4 

0    3 

2 

0    3 

0 

0    3 

0 

0    3 

1 

0    3 

3 

0    3 

3 

0    3 

'  6 

0.  4 

15 

0    4 

20107 

82    9 

121 

83    4 

7 

83    5 

6 

83    5 

1 

83    5 

0 

83    5 

0 

83    5 

0 

83    5 

0 

83    5 

0 

83    5 

25 

83    6 

1 

83    6 

3 

83    6 

1 

83    6 

0 

83    6 

0 

83    6 

0 

83    6 

3 

83    6 

2 

83    6 

9 

83    7 

1027 

87    9 

5 

87    9 

1 

87    9 

0 

87    9 

3 

87    9 

1 

87    9 

0 

87    9 

0 

87    9 

0 

87    9 

1 

87    9 

12 

88    0 

2929 

100    0 

r,g 


Table  XIV 
Global  Magnetic  Cutoff  Values 


L-Value 


Bcutoff 


1.0949 

[J.23901 

1.1152 

22998 

1.1350 

22182 

1.1559 

22288 

1.1752 

22353 

1.2628 

24017 

1.4817 

26770 

1.8331 

27573 

2.3168 

.29762 

2.9330 

.33311 

3.6816 

37001 

4.5626 

.40274 

5.5760 

.43178 

6.7218 

.45502 

8.0000 

.46951 

Note:  The  Table  may  be  interpolated  for  other  L-ral-aes  (Ref.  16). 
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The  overall  flow  for  the  program  is  as  follo-ws  (see  Appendix  E): 

(1)  Use  an  initial  guess  for  the  fundamental  decay  time  and  the  theoretical  ratios 

of  the  higher-mode  dc~y  times  to  compute  the  starting  decay  tin  es.  The  first 
guess  at  fundamental  decay  time  is  shown  in  Table  XV  as  a  function  of  L  and 
energy.  This  table  is  interpolated  in  L  and  E  for  starting  value. 

(2)  Numerically  compute  the  integrals  of  the  functions  shown  in  Appendjx  C  at  oach 

data  point  using  the  initial  guessed  decay  times.  The  integration  s  performed 
by  a  12-point  Gauss/Legcndre  Quadrature  at  each  data  point  for  e:ich  function. 
Only  the  fundamental  decay  mode  and  two  higher  modes  were  c(  mputod. 

(3)  Compute  the  linear  amplitude  constants  using  a  least-squares  matrix-solving 
routine.  Compute  the  standard  error  of  the  flt.  The  individual  points  are 
weighted  by  the  reciprocal  of  their  fractional  errors,  which  have  been  normal- 
ised to  an  average  value  of  one. 

(4)  Optimize  the  initially  computed  yeut  value  by  using  a  "Golden  Section  Min- 

imim"  function.  This  section  attempts  to  let  the  data  determine  the  "best" 
cutoff  value.  If  the  data  are  not  distributed  down  toward  the  loss  cone,  the 
Golden  Minimum  function  will  provide  an  erroneous  ycul  value  Hence,  the 
function  specifies  that  the  computed  yeut  must  be  within  ±10  percent  of  the 
initial  input  value.  If,  in  the  iteration  process,  the  function  tries  to  exceed  these 
limits,  the  initial  guess  of  yeut  is  specified  as  the  correct  value. 

(5)  Repeat  the  *leps  (2)  and  (3)  with  the  optimum  yeut  value. 

(6)  To  determine  the  best  combination  of  decay  times  and  amplitude  constants, 
the  decay  times  are  changed  by  a  "random-walk"  method,  either  increased 
or  decreased.  The  amount  by  which  the  decay  times  may  randomly  vary  is 
determined  by  using  a  cumulative  normal  probability  distribution  about  the 
initial  value. 

The  standard  deviation  used  was  estimated  according  to  this  author's  best  guess  of 
the  accuracy  of  the  tables  and  curves  shown  by  Stassinopoulos  (Ref.  53:31-32,  40-44)  and 
West  (Ref.  63:50-54).  The  Stassinopoulos  tables  do  not  always  correlate  with  each  other 
or  with  the  Stassinopoulos  curves,  and  the  West  curves  do  not  always  correl  ite  with  the 
Stassinopoulos  in  the  L-regions  of  overlap.  The  initial  fundamental  decay  times  in  Table  XV 
were  computed  by  this  author  by  interpolation  of  the  tables  and  curves  in  Reference  53  and 
by  reading  and  extrapolation  of  the  curves  in  Reference  83.  The  author's  estimates  of  their 
standard  deviations  are  shown  in  Table  XVI.  Of  particular  note  is  the  very  lov  confidence 
in  decay  times  at  high  L  values  (above  2.2)  and  at  high  energies  (above  2.0  ivleV).  There 
have  been  few  computations  of  decay  times  in  these  regions,  and  the  results  vary  widely 
(Rcfs.  53,  03). 

The  effect  of  using  the  cumulative  normal  probability  distribution  and  the  above  stan- 
dard deviations  is  to  restrict  the  "random-walk"  of  the  decay  times  away  from  their  initial 
values.  The  farther  away  from  the  initial  "mean  value"  that  the  random  val  le  generator 
places  the  new  decay  time,  the  greater  the  "push"  that  is  created  back  toward  he  mean  for 
the  next  walk.  The  distance  away,  however,  is  defined  in  terms  of  the  standa.d  deviation, 
so  that  r  values  in  which  there  is  low  confidence  may  vary  by  larger  amounts  than  may  r 
values  in  which  there  is  greater  confidence.  The  r's  are  also  restricted  by  the  co  istraint  that 
ro  >  r\  >  r2. 
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With  each  change  of  decay  times,  steps  (2)  and  (3)  are  repeated  up  to  90  times,  in  an 
attempt  to  get  a  "better"  Qt  by  reducing  the  standard  error  of  the  fit. 
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Table  XVT 

Percent  Standard  Deviation  of  Initial 
Fundamental  Decay  Times  (in  Table  XV) 


lr  Range 

1.02- 

1.21- 

1.1 1— 

1.81- 

> 

1.20 

1.40 

1.80 

2.20 

2.20 

Energy  Range 

MeV 

E   <   0.0 

30 

30 

20 

35 

200 

0.0  <  E   <   2.0 

30 

20 

20 

75 

100 

E   >  2.0 

50 

20 

50 

100 

100 

(estimated  from  Ref.  53:31-32,  40-44  and  R.?f.  53:50-54) 


44 


(7)  If  a  better  fit  is  found  in  less  than  98  random  walks,  and  if  less  than  32  "better" 
values  hare  been  found,  steps  (2),  (3)  and  (0)  are  repeated  up  to  32  times,  using 
the  latest  "best"  decay  times. 

(8)  Up  to  five  repeats  of  steps  (4)  and  (5)  may  be  carried  out,  but  only  if  32  "better" 
fits  are  found.  Otherwise,  only  the  first  yeut  optimization  is  performed. 

(9)  The  program  will  always  end  at  the  point  where  96  random-walk3  have  failed 
to  reduce  the  error  of  the  fit  from  the  previous  "best-fit".  The  iccay  times, 
yeut,  ^nd  amplitude  constants  are  printed  for  each  "better"  fit,  with  the  final 
values  being  those  of  the  "best"  fit. 

(10)  After  the  best  lit  is  found,  the  program  uses  the  fitted  constants  and  decay 
times  to  compute  fluxes,  at  specified  times,  from  the  equator  to  the  best  cutoir 
value.  These  specified  times  correspond  to  default  times  plus  1/2  ime  window 
in  ^he  program  which  plots  raw  data.  Thus,  the  computed  (luxe:-  are  used  to 
generate  curves  which  are  plotted  over  raw  dat3  points  at  approximately  the 
same  times. 

Variations  in  the  program  were  tried  in  order  to  study  the  effects  of  the  fitting  process. 
The  principal  modifications  were  in  the  random  decay-time  variations. 

The  first  modification  was  to  fix  the  ratios  of  the  r's  to  their  theoretical  values  of 
1:1/9:1/25  (Ref.  48:16-19).  The  effect  of  this  modification  is  to  significantly  reduce  the 
degrees  of  freedom  available  to  fit  the  data. 

The  second  modification  was  a  compromise  between  the  above  and  complete  freedom  to 
random  walk.  This  variation  allowed  r0  and  r^  to  vary  freely  (subject  to  the  usual  contraint 
that  ro  >  n),  but  T\  and  u  were  fixed  at  their  theoretical  ratio  of  25:9.  The  principal  reason 
for  this  attempt  was  the  very  short  time-span  covered  by  some  satellite/burst  combinations. 
For  example,  the  Russian  2  Burst  covered  only  4  1/2  days  prior  to  the  Russian  3  injection. 
If  the  data  cover  a  time-span  which  is  small  compared  to  the  fundamental  decay  constant, 
then  the  Gtted  value  of  that  decay  time  may  be  suspect.  However,  if  the  time-s oan  is  not  so 
small  compared  to  the  decay  times  of  the  higher  modes,  one  should  have  more  confidence 
in  the  fitted  values  of  those  decay  time.  Stated  another  way,  the  fit  may  not  b<-  sensitive  to 
the  fundamental,  but  may  still  be  sensitive  to  the  higher  modes. 

In  this  analysis,  not  all  L-values  covered  by  the  satellite  in  Table  II  were  examined. 
Over  1000  plots  of  flux  versus  x  were  examined  at  various  L-values  for  the  satellite/burst 
combinations  of  Tables  I  and  II.  In  many  cases,  the  satellite  orbits  were  of  such  high 
inclination  (toward  polar)  that  data  were  not  available  close  to  equatorial  B- values.  Because 
of  this,  and  because  of  the  very  large  volume  of  data,  those  satellite/burst/L  combinations 
which  provided  the  most  consistent  coverage  from  the  equator  to  the  loss  cone  were  chosen 
for  study. 

Generally,  the  data  fell  into  three  regions  of  x-spaee  (with  some  overlap),  legion  I  was 
from  the  equator  to  the  point  where  flux  dropped  off  sharply  into  the  loss  cone  (x  =  0  to 
i  rs  0.6  or  0.8).  Region  II  was  the  area  where  the  flux  "turned  the  corner"  from  a  rather 
"horizontal"  curve  to  a  rather  "vortical"  curve  into  the  loss  cone  (x  rs  0.7  to  0.'.1).  Region  Ml 
was  the  loss  cone  region,  where  the  flux  levels  dropped  sharply  to  the  cutoff  v  lue  [x  sa  .9 
for  a  typical  L- value). 


" 


Representative  plots  of  raw  data  in  the  three  regions  are  shown  in  Figures  4  through  8. 
These  plots  are  typical  in  that  data  seldom  cover  all  tbree  regions  in  a  short  tinrvspan.  Gaps 
in  the  data  appear  because  of  satellite  orbital  coverage  and  because  of  lack  of  on-satellite 
recorders.  When  a  satellite  lacked  a  recorder,  data  were  collected  on'.7  as  the  satellite  passed 
above  the  radio  horizon  of  a  ground  tracking  station.  Figure  Q  is  an  eight-hour  plot  of  typical 
satellite  coverage  over  all  B-L  space. 
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0.5 


X  =  (1-Beq/B)1/2 

Figure  4(a).  Raw  Allotxette  Data  Collected  Only  in  Region  ITT  (at  buret  tine). 
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BURST  :  RUSSIAN  2  SAT  :  ALL0UETTE1 

TIME  :     3  (dy)      0  (hr)      0.00  (min)  3.000(days) 

L  :  2.50      Beq  :  .020      BCUT  :  0.515      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  4(b).  Raw  Allonette  Data  Collected  Only  in  Region  HI  (3  days  post-burst). 
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BURST  :  STARFSH  SAT  :  TELS~AR 

TIME  :     1  (dy)      0  (hr)      0.00  (min)  1.00C(doys) 

L  :  1.80      3eq  :  .053      BCUT  :  0.470      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gcuss)  :  1.000 
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LEGEND 
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X  =  (l-3eq/B)V2 

Figure  5.  Ra-w  Telstar  Data  Collected  in  Regions  I  and  IH  (early  time). 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :    15  (dy)      0  (hr)      0.00  (min)  15.000(days) 

L  :  1.70      Beq  :  .063      BCUT  :  0.460      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 


X  =  (1-Beq/B)1/2 

Figure  6.  Raw  Telstar  Data  Collected  in  Region  III  (mid- range  time). 
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Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  7.  Raw  Telstar  Data  Collected  in  Regions  1,  II,  and  HI  (mid-range  Ume) 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :  100  (dy)      0  (hr)      0.00  (min)         100.000(days) 
L  :  1.80      Beq  :  .053      BCUT  :  0.470      TW  (days)  :     5.0 
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Figure  8(a).  Raw  Telstar  Data  in  Regions  I  and  TI  (late  time). 
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Figure  8(b).  Raw  Telstar  Data  in  Regions  I,  II,  and  III  (late  time). 
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AFWL     TRAPPED    ELECTRON     DATA       BASE 
BURST  :  STARFISH  SATELLITE  :  TELSTAR 

DAY  :  237    HOUR  :    1    MINI  :    0  #  :     98 

CHANNEL  NUMBER  :  4      TIME  WINDOW  :    8.00 
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Figure  9.  Typical  Raw  Data  Coverage  for  8  Hours  [B-L  values  not  shown). 


8 


r>4 


Figures  10  (a)  through  (g)  shows  raw  data  from  Explorer  XV  over  the  time  period  of  the 
Russian  2  and  3  bursts,  at  times  corresponding  to  those  shown  in  Figure  3,  for  comparison. 
It  is  apparent  that  fewer  points  are  plotted  from  the  AFWL  data  base  in  Figures  10  (a) 
through  (g)  than  were  plotted  by  Roberts  in  Figure  3  for  the  1.9  Mc~  electrons,  even  though 
the  satellite  and  time  periods  of  coverage  were  the  same.  These  discrepancies  arc  not  resolved 
in  this  study,  although  the  previously  alluded-to  L-rounding  procedure  may  account  for  the 
differences. 

The  temporal  progression  of  flux  after  one  burst  from  one  satellite  (Telstar),  at  a  typical 
L- value  of  1.9,  is  shown  in  Figures  11  (a)  through  (h).  Again,  typical  gaps  in  coverage  of 
Regions  I,  II,  and  III  are  readily  apparent  over  the  90  days  shown. 

As  previously  noted,  only  fundamental  and  two  higher  eignemodes  (/o,/i,/2)  are  used 
to  Ot  the  satellite  data  (Appendix  C).  The  third  eigenmode  (/:»)  is  not  used  becai  se  Shulz  and 
Lanzerotti  (Ref.  45:163)  suggest  that  the  modes  greater  than  two  rapidly  vanish.  Initially,  the 
steady-state  function  (/oo)  was  used  along  with  the  fundamental,  but  this  solution  was  found 
to  be  "competing"  with  the  fit  for  the  fundamental,  in  the  sense  that  the  linear  constants 
tended  to  be  roughly  equal  in  magnitude  and  tended  to  alternate  in  sign.  Shulz  (Rcf.  47:22- 
23)  indicates  that  the  /oo  and  /o  functions  should  resemble  each  other  in  shape.  Also,  PQtser 
(Ref.37)  has  subtracted  the  background  from  all  flux  data,  and  the  background  should 
approximate  the  steady-state,  /oo,  solution.  For  these  reasons,  the  steady-state  solution  was 
not  used  in  the  data  fit. 
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BURST  :  RUSSIAN  2  SAT  :  EXPLORER15/EXPL0RER15 

TIME  :     0  (dy)      4  (hr)     58.07  (min)  0.207(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  10(a).  Raw  Explorer  15  Data  at  Times  Corresponding  to  Figure 
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BURST  ;  RUSSIAN  2  SAT  :  EXPI0RER15/EXPL0RER15 

TIME  :     2  (dy)     19  (hr)     22.09  (min)  2.807(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     2.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  10(b).  Ra-w  Explorer  15  Data  at  Times  Corresponding  to  Figure  3. 
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BURST  :  RUSSIAN  2  SAT  :  EXPLORER15/EXPLORER15 

TIME  :     9  (dy)     19  (hr)     22.09  (min)  9.807(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     2.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  10(c).  Raw  Explorer  15  Data  at  Times  Corresponding  to  Figure 
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BURST  :  RUSSIAN  2  SAT  :  EXPL0RER15/EXPL0RER15 

TIME  :    20  (dy)     19  (hr)     22.09  (min)  20.807(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     2.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  lOfd).  Raw  Explorer  15  Data  at  Times  Corresponding  to  Fignre  3. 
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BURST  :  RUSSIAN  2  SAT  :  EXPLORER15/EXPIORER15 

TIME  :    29  (dy)     19  (hr)     22.09  (min)  29.807(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     4.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  10(e).  Raw  Explorer  15  Data  at  Times  Corresponding  to  Figure  3. 
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Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  10(f).  Raw  Explorer  15  Data  at  Times  Corresponding  to  Figure  3. 
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BURST  :  RUSSIAN  2  SAT  :  EXPL0RER15/EXPL0RER15 

TIME  :    63  (dy)     19  (hr)     22.06  (min)  63.807(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :    10.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  10(g).  Raw  Explorer  15  Data  at  Times  Corresponding  to  Figure  3. 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :     1  (dy)      0  (hr)      0.00  (min)  1.000(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  11(a).  Temporal  Progression  for  Ra-w  Telstar  Data  from  1  to  90  Days  Post-Burst. 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :     2  (dy)      0  (hr)      0.00  (min)  2.000(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  11(b).  Temporal  Progression  for  Raw  Telstar  Data  from  1  to  90  Dap-,  Post-Burst. 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :     4  (dy)      0  (hr)      0.00  (min)  4.000(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  11(c).  Temporal  Progression  for  Raw  Telstar  Data  from  1  to  90  Day;  Post-Burst. 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :     6  (dy)      0  (hr)      0.00  (min)  6.000(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     2.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  11(d).  Temporal  Progression  for  Raw  Telstar  Data  from  1  to  00  Days  Post-Burst. 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :    10  (dy)      0  (hr)      0.00  (min)  10.000(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  11(e).  Temporal  Progression  for  Raw  Telstar  Data  from  I  to  90  Day.-,  Post- Burst. 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :    15  (dy)      0  (hr)      0.00  (min)  15.000(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  11(f)-  Temporal  Progression  for  Rn  Telstar  Data  from  1  to  90  Days  Post-Burst. 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :    60  (dy)      0  (hr)      0.00  (min)  60.000(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  11(g).  Temporal  Profession  for  Raw  Telstar  Data  from  1  to  90  Days  Post-Burst. 
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BURST  :  STARFISH  SAT  :  TELSTAR 

TIME  :   90  (dy)      0  (hr)      0.00  (min)  90.000(days) 

L  :  1.90      Beq  :  .045      BCUT  :  0.479      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  11(h).  Temporal  Progression  for  Raw  Telstar  Data  from  1  to  90  Day:;  Post-Burst. 
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V.  RESULTS  AND  DISCUSSION 

Figures  12  and  13  show  the  typical  temporal  progression  of  raw  data  at  two  representative 
L  values,  L  =  2.3  and  L  =  2.1.  Included  on  those  figures  are  the  [lux  curves  calculated  from 
the  amplitude  coefficients  and  decay  times  fitted  to  that  data  by  Electrofit.  As  previously 
noted,  the  curves  shown  on  any  given  figure  may  vary  slightly  in  time  from  the  individual 
data  points,  since  the  curves  represent  only  the  one  time-instant  corresponding  to  the  start 
time  plus  one-half  the  time  window  of  the  plot.  Figures  12  and  13  represent  maximum 
freedom  of  the  fit,  since  each  r  is  free  to  vary  independently  of  the  other  r's. 

Figure  14  shows  the  same  temporal  data  progression  with  different  fitted  curves.  The  curves 
in  this  figure  were  calculated  under  the  restriction  that  the  ratios  of  the  r's  ^vere  fixed  at 
the  theoretical  values  of  1:1/9:1/25;  hence,  only  T\  was  really  varying  in  a  random  sence. 

Figure  15  shows  similar  data  progression  and  fitted  curves  where  the  r's  wer<>  allowed  to 
freely  vary,  subject  to  the  constraint  that  tq  >  rx  >  r2.  Also,  this  was  an  initial  -un  with  500 
tries  for  convergence,  vice  98  tries.  In  addition,  the  computed  flux  was  used  in  the  iteration 
process  to  eliminate  non-phys'cal  (negative)  solutions.  The  solution  was  restricted  to  being 
positive  at  four  different  times  and  at  five  different  equatorial  pitch  angles.  Alt  lough  this  is 
an  initial  run  with  only  a  small  amount  of  data  used,  it  indicates  that  the  recommendations 
in  Chapter  VI  should  be  pursued. 

Table  XVII  shows  the  fitted  linear  amplitude  coefficients  and  exponential  decay  times  (in 
days)  at  constant  Energy  and  //-shell  for  a  representative  subset  of  the  data  examined  in 
this  study.  The  sheer  volume  of  resultant  fits  and  calculated  curves  made  it  impossible  to 
include  all  results  in  this  report.  The  values  shown  in  Table  XVII  were  computed  with  the 
r's  allowed  to  vary  freely. 

Figures  12  and  13  are  considered  representative  of  the  results  of  the  fitting  model.  It  is 
apparent  from  examination  of  Figures  12  and  13  that  the  pitch-angle  diffusion  theory 
and  flux  computation  method  developed  in  this  report  are  consistent  with  experimental 
measurements.  Higher  eigenmodes  are  readily  discernible  at  the  lowest  energy  in  the  Figures 
or  the  first  five  days,  and  are  not  visible  at  later  times.  These  higher  eigeumodes  are  not 
discernible  in  the  plots  at  the  higher  energies;  Table  XVII  shows  that  theso  modes  are 
present,  but  at  much  lower  levels,  in  some  cases.  The  calculated  curves  show  a  reasonable 
fit  to  the  data  even  at  late  times,  consistent  with  the  assumption  of  exponenti  il  decay. 

Examination  of  Figures  14  and  15  shows  that  some  inconsistencies  remain  to  bo  resolved  in 
the  use  of  this  fitting  method.  The  caluclated  flux  curves  for  the  low-energy  data  points  also 
appear  to  be  slightly  lower  than  optimum,  while  the  high-energy  data  appear?  to  be  fitted 
better  by  the  model.  This  may  be  partially  explained  by  the  existence  of  many  lower  data 
points  in  the  vicinity  of  the  loss  cone  which  tend  to  "drag  down"  the  overall  fit  Any  change 
in  the  cutoff  value  significantly  affects  the  fit  by  shifting  the  "vertical"  portion  of  the  curve 
toward  or  away  from  those  points.  However,  this  cannot  account  for  the  entire  difference 
of  the  curves  from  the  data.  The  energy  dependence  may  be  treated  less  tlmn  optimally 
in  the  choice  of  functional  solution  made  in  this  study.  If  the  chosen  functional  solution  is 
not  the  correct  one,  it  may  be  nevertheless  close  enough  to  correct  to  give  a  reasonable  fit 
at  high  energy,  but  be  significantly  in  error  at  lower  energy.  A  different  functional  solution 
may  fit  the  lower-energy  groups  better  and  yet  still  fit  the  high-energy  elect  ons  as  well. 
Recommendation  for  the  use  of  other  functions  is  made  in  Chapter  VI. 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     0  (dy)      0  (hr)      0.00  (min)  0.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     0.5 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  12(a).  Temporal  Profession  of  Raw  Data  and  Fitted  Flux  Curves  at  L  —  2.3  for 
Russian  3  Burst  (rn  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     0  (dy)     12  (hr)      0.00  (min)  0.500(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     0.5 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  12(b).  Temporal  Progression  of  Raw  Data  3nd  Fitted  Flux  Curves  at  L  =  2.3  for 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     1  (dy)       0  (hr)      0.00  (min)  1.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  12(c).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.3  for 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     2  (dy)      0  (hr)      0.00  (min)  2.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  12(d).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L 
Russian  3  Bnrst  [r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     3  (dy)      0  (hr)      0.00  (min)  3.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  12(e).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.3  for 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     4  (dy)      0  (hr)       0.00  (min)  4.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  12(f)-  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     5  (dy)      0  (hr)      0.00  (min)  5.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  12(g).  Temporal  Progression  of  Ra-w  Data  and  Fitted  Flux  Curves  at  /,  =  2.3  for 
Russian  3  Burst  (rn  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :    25  (dy)      0  (hr)      0.00  (min)  25.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  12(h).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.3  for 
Russian  3  Burst  (rrt  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :    60  (dy)      0  (hr)      0.00  (min)  60.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Fignre  12(1)-  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curres  at  L 
Russian  3  Burst  (r„  freely  varying). 


=  2.3  for 


10 


,io 


10" 


10° 


--    10' 


10" 


10' 


10' 


sn 


10'°  p 


(J 

Ld 
00 


(J 

o 

o: 

h- 
(J 
Ld 
_J 
LJ 

X 

3 


o 

o: 
k- 
O 

LJ 
_J 
Ld 

_J 
< 

O 
Ld 


1CT 


BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME:     0  (dy)      0  (hr)      0.00  (min)  O.OOO(days) 

L  :  2.40      Beq  :  .022      BCUT  :  0.511      TW  (days)  :     0.5 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  13(a).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2A  for 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     0  (dy)     12  (hr)       0.00  (min)  0.500(days) 

L  :  2.40      Beq  :  .022      BCUT  :  0.511      TW  (days)  :     0.5 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  13(b).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.4  for 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [0-15] 

TIME  :     1  (dy)      0  (hr)      0.00  (min)  1.000(days) 

L  :  2.40      Beq  :  .022      BCUT  :  0.511      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  13(c).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  —  2.4  for 
Russian  3  Burst  (rn  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     2  (dy)      0  (hr)      0.00  (min)  2.000(days) 

L  :  2.40      Beq  :  .022      BCUT  :  0.511      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  13(d).  Temporal  Progression  of  Raw  Data  and  Fitted  Flax  Curves  at  L  =  2.4  for 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     3  (dy)      0  (hr)      0.00  (min)  3.000(days) 

L  :  2.40      Beq  :  .022      BCUT  :  0.511      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  13(e).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  -  2.4  for 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :     4  (dy)      0  (hr)      0.00  (min)  4.000(days) 

L  :  2.40      Beq  :  .022      BCUT  :  0.511      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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FifT-irc  13(f).  Temporal  Progression  of  Raw  Data  and  Fitted  Flax  Curves  at  L 
Russian  3  Burst  (r„  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :    25  (dy)      0  (hr)      0.00  (min)  25.000(days) 

L  :  2.40      Beq  :  .022      BCUT  :  0.511      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  13(g).  Temporal  Progression  of  Rrx  Data  and  Fitted  Flux  Curves  at  I  =  2.1  Tor 
Russian  3  Burst  (r„  froely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 

TIME  :    60  (dy)      0  (hr)      0.00  (min)  60.000(days) 

L  :  2.40      Beq  :  .022      BCUT  :  0.511      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  13(h).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L 
Russian  3  Burst  (rn  freely  varying). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [052] 

TIME  :     0  (dy)      0  (hr)      0.00  (min)  O.OOO(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     0.5 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  14(a).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Currcs  at  L  =  2.3  for 
Russian  3  Burst  (r„  fixed  1:1/9:1/25). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [052] 

TIME  :     0  (dy)     12  (hr)      0.00  (min)  0.500(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     0.5 
Bmln  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  14(b).  Temporal  Progression  of  Row  Data  and  Fitted  Flrnc  Carres  at  /,  =  2.3  for 
Russian  3  Borst  (r„  toed  1:1/0:1/25). 
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BURST  :  RUSSIAN  3  SAT  :  EXPLORER15  [052] 

TIME  :     1  (dy)      0  (hr)      0.00  (min)  LOOO(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  14(c).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Cuttcs  at  L  =  2.3  for 
Russian  3  Burst  (r„  fixed  1:1/9:1/25). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [052] 

TIME  :     2  (dy)      0  (hr)      0.00  (mm)  2.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  14(d).  Temporal  Profession  of  Raw  Data  and  Fitted  Flux  Curves  at  h  —  2.3  for 
Russian  3  Burst  (r„  fixed  1:1/9:1/25). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [052] 

TIME  :     3  (dy)      0  (hr)      0.00  (min)  3.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  14(e).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.3  for 
Russian  3  Burst  (r„  fixed  1:1/9:1/25). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [052] 

TIME  :     4  (dy)      0  (hr)      0.00  (min)  4.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000     • 
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Figure  14(f).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curres  at  L 
Russian  3  Burst  (rn  fixed  1:1/0:1/25). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [052] 

TIME  :    50  (dy)      0  (hr)      0.00  (min)  50.000(days) 

L  :  2.30      Beq  :  .025      BCUT  :  0.506      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 


LEGEND 
a  =  .500  MEV 
a  =  .500  MEV 


10" 

0  0.5 

X  =  (1-Beq/B)1/2 

Figure  14(g).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.3  for 
Russian  3  Burst  (rn  fixed  1:1/9:1/25). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :     0  (dy)      0  (hr)      0.00  (min)  0.000(days) 

L  :  2.40       Beq  :  .022       Bcut  :  0.302       TW  (days)  :     0.5 

Bmin  (Gauss)  :  0.001        Bmax  (Gauss)  :  1.000 
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Figure  15(a).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.1  for 
Russian  3  Burst  (rra  freely  varying  with  phys-cal  oontraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :     0  (dy)     12  (hr)      0.00  (min)  0.500(day3) 

L  :  2.40       Beq  :  .022       Bcut  :  0.302       TW  (days)  :     0.5 
Bmin  (Gauss)  :  0.001        Bmax  (Gauss)  :  1.000 
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Figure  15(b).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L 
Russian  3  Uurst  (r„  freely  varying  with  physical  contraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :     1  (dy)      0  (hr)      0.00  (min)  1.000(days) 

L  :  2.40       Beq  :  .022      Bcut  :  0.302      TW  (days)  :     1.0 

Bmin  (Gauss)  :  0.001        Bmax  (Gauss)  :  1.000 
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Figure  15(c).  TempornI  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.4  for 
Russian  3  Burst  (rn  freely  varying  with  physical  c.ontraints). 
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BURST  •  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :     2  (dy)      0  (hr)       0.00  (min)  2.000(days) 

L  :  2.40       Beq  :  .022       Bcut  :  0.302      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001        Bmax  (Gauss)  :  1.000 
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Figure-  15(d).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Carves  at  L 
Russian  3  Durst  (r„  frcciy  varying  with  physical  contraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :     3  (dy)      0  (hr)      0.00  (min)  3.000(days) 

L  :  2.40      Beq  :  .022      Bcut  :  0.302      TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  15(e).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Carres  at,  L 
Russian  3  Burst  (r„  freely  varying  wUa  physical  coatraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPLORERS  [568  j 

TIME  :     4  (dy)      0  (hr)      0.00  (min)  4.000(doys) 

L  :  2.40       Beq  :  .022       Bcut  :  0.302       TW  (days)  :     1.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  15(f).  Temporal  Progression  of  Ra-w  Data  and  Fitted  Flux  Curves  at  L 
Russian  3  Burst  (r„  freely  varying  with  physical  cootraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :     5  (dy)      0  (hr)      0.00  (min)  5.000(days) 

L  :  2.40      Beq  :  .022       Bcut  :  0.302       TW  (days)  :     1.0 

Bmin  (Gauss)  :  0.001        Bmax  (Gauss)  :  1.000 
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Figure  15(c).  Temporal  Progression  of  Raw  Data  and  Fitted  Flnx  Carres  at  L  =  2.4  for 
Russian  3  Burst  (rn  freely  varying  with  physical  eonlraiuts). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :     6  (dy)      0  (hr)      0.00  (min)  6.000(days) 

L  :  2.40      Beq  :  .022       Bcut  :  0.302      TW  (days)  :     2.0 
Bmin  (Gauss)  :  0.001        Bmax  (Gauss)  :  1.000 
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Figure  15(h).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  —  1A  for 
Russian  ?,  Burst  (rn  freely  varying  with  physical  contraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :     8  (dy)       0  (hr)      0.00  (min)  8.000(days) 

L  :  2.40      Beq  :  .022       Bcut  :  0.302       TW  (days)  :     2.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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FiRTire  15(i).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.1  for 
Russian  ?>  Bxirst  (rn  freely  varying  with  physical  contraints). 


tn-i 


10 


10 


u 

LxJ 


(J 

S 

00 

o 

en 

t— 

(J 

Ld 


X 

D 


Z 

o 

i— 
o 


UJ 

_i 
< 

ct: 
o 

UJ 


10' 


10° 


10' 


Or 


10* 


10= 


itf 


BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :    10  (dy)      0  (hr)       0.00  (min)  10.000(days) 

L  :  2.40      Beq  :  .022       Bcut  :  0.302      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001        Bmax  (Gauss)  :  1.000 
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Figure  15(j).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =  2.4  for 
Russian  o  Burst  (r„  freely  varying  with  physical  contraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :    35  (dy)      0  (hr)      0.00  (min)  35.000(days) 

L  :  2.40      Beq  :  .022      Bcut  :  0.302      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  15(k).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  3t  //  =  1A  for 
Russian  3  Burst  (r„  freely  varying  -with  physical  contraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :    50  (dy)      0  (hr)      0.00  (min)  50.000(days) 

L  :  2.40       Beq  :  .022       Bcut  :  0.302      TW  (days)  :     5.0 
Btnin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  15(1).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  //  =  2.4  for 
Russian  3  Burst  (rn  frocly  varying  with  physical  contraints). 
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BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [568] 

TIME  :    60  (dy)      0  (hr)      0.00  (min)  60.000(doys) 

L  :  2.40      Beq  :  .022      Bcut  :  0.302      TW  (days)  :     5.0 
Bmin  (Gauss)  :  0.001       Bmax  (Gauss)  :  1.000 
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Figure  15(m).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  rst  L  =  2.4  Tor 
Russian  i  Burst  (rn  freely  vnryinp  with  physical  contraints). 
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The  short  time-spans  covered  by  most  of  the  data  sets  in  the  Trapped  Electron  Data  Base 
have  not  allowed  adequate  determination  of  n>.  The  values  of  t\)  presented  in  Table  XVII 
arc  suspect  for  this  reason.  However,  more  confidence  is  held  in  the  values  >f  rt  and  r2 
computed  by  this  study.  In  most  data  sets  studied,  the  time  period  of  the  d.tta  was  long 
enough  compared  to  the  values  of  rt  and  r^  to  compute  reasonable  values.  More  consistency 
between  data  sets  was  found  for  r,  and  r2  than  was  found  for  r0.  Of  course,  then  rn  may  be 
simply  computed  by  multiplying  r,  by  the  theoretical  ratio  corresponding  lo  the  solution. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

The  pitch-angle  diffusion  theory  and  flux  calculation  method  developed  in  this  study 
are  reasonably  consistent  with  oxperimental  data  in  most  cases.  The  method  si  ows  promise 
as  a  way  of  calculating  improved  decay  times  over  a  broad  region  of  the  ma 'netosphere, 
particularly  for  higher  energy  electrons.  Additionally,  a  comprehensive  summary  table  of 
the  best  available  literature  values  for  ro  has  been  prepared  (Table  XV).  A  complete  table 
covering  such  a  broad  range  of  energy  and  L  has  not  heretofore  been  compiled. 

Time  and  computer  limitations  have  precluded  examination  of  the  entire  data  base. 
However,  the  overall  efficiency  of  the  model  developed  in  this  study  has  been  demonstrated.  It 
is  expected  that  improved  predictions  of  fluxes  and  improved  inputs  to  the  AFW,  SPECTER 
Codes  (Ref.  7-9)  will  be  available  as  a  result  of  this  work.  This  will  enhance  tie  capability 
for  calculation  of  satellite  operational  environments,  and  will  improve  the  abili  y  to  predict 
satellite  survivability  and  vulnerability. 


Recommendations 

The  following  recommendations  for  further  study  are  proposed  to  improve  the  model 
developed  in  this  study,  and  to  resolve  the  problem  of  non-physical  solutions  for  some  data 
sets.  At  the  direction  of  this  author,  several  of  these  recommendations  are  presently  being 
implemented  at  the  Air  Force  Weapons  Laboratory. 

(1)  Incorporate  an  estimation  of  the  error  in  the  magnetic  field  value,  D,  of 
individual  data  points.  This  corresponds  to  a  "horizontal"  error  bar  on  the 
plots  such  as  Figure  12.  It  is  known  that  errors  exist  in  the  compulation  of  the 
B  value  associated  with  each  data  point,  and  this  error  would  prove  significant 
in  the  region  of  the  loss  cone.  Improved  weighting  of  the  data  points  with  the 
combined  errors  should  improve  the  flt. 

(2)  Utilize  the  //-rounding  procedure  of  Chapter  IV  to  add  data  points  to  each  L- 
value  considered.  This  should  aid  in  filling  gaps  in  the  data,  whic'i  necessarily 
will  improve  the  fit. 

(3)  Discard  the  highest  eigenraode  if  the  fitted  value  of  its  decay  tim'>  is  less  that 
0.25  day,  and  fit  the  data  with  only  the  fundamental  and  one  higher  mode. 
A  value  of  0.25  day  is  not  reasonably  fitted  within  the  errors  in  the  Trapped 
Electron  Data  Base.  This  mode  would  be  essentially  lost  within  th;  first  day  in 
any  case. 

(4)  Perform  an  "energy-scaling"  of  the  data  to  combine  data  from  different  satel- 
lites. This  involves  assuming  some  form  of  energy  spectrum,  such  as  e~iE^l''\ 
including  an  estimate  of  the  error  in  the  assumed  value  of  Ko-  The  new  flux 
would  be  simply  the  old  flux  multiplied  by  a  factor  e~[Bn*'~E'u),E,k .  If  the 
error  in  the  old  flux  were  50  percent  (a  value  consistent  with  the  Data  Rase) 
and  the  error  in  Eq  were  25  percent,  and  if  £o  ^  1-1  (a  fission  spectrum),  then 
the  resultant  error  in  the  shifted  flux  would  be  less  than  1  percent  different 
from  that  of  the  old  flux,  when  shifting  from  660  KeV  to  500  KeV.  Thus,  it 
should  be  possible  to  combine  data  sets  and  fill  the  gaps  in  data  to  achieve  a 
better  fit.  However,  the  limitation  of  the  assumed  spectrum  still  exists. 
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(5)  Utilize  only  late-time  data  where  it  is  available,  3nd  Gt  only  the  fundamental 
mode  to  it,  under  the  assumption  that  higher  modes  will  have  decayed  away. 
Then,  with  the  fundamental  mode  constrained,  go  back  to  early  time  data  aod 
attempt,  to  fit  higher  modes,  if  any.  This  recommendation  is  of  limited  value 
because  only  two  data  sets  cover  a  significant  period  of  time:  Starfish/Telstar 
and  Russian  3/Explorer  5.  However,  it  may  provide  a  better  value  of  r0  than  is 
presently  available  at  some  h  values. 

(ft)  Revise  the  Gtting  functions  to  other  Bessel  function  solutions  thnn  the  J-\ji 
solution  and  compare  results  This  was  performed  only  on  one  oth<*r  functional 
solution,  Jo,  with  only  one  small  data  set,  and  the  results  showed  no  improve- 
ment over  the  present  study.  However,  to  be  certain  that  the  an  hor  has  not 
simply  chosen  the  wrong  functional  form,  more  functional  forms  must  be  in- 
vestigated. 

(7)  In  the  fitting  program,  constrain  the  directional  flux  to  physical  reality  at  each 
iteration  (i.e.,  compute  flux  at  each  iteration  and  reject  non-physical  solutions 
at  each  iteration,  rather  than  after  converging  on  a  final  solution).  This  should 
converge  to  a  physically  real  solution  if  it  converges  at  all;  however,  significant 
increases  in  computer  processing  time  may  accrue. 

(8)  Increase  the  number  of  tries  for  a  solution  from  98  to  200  and  reduce  the  number 
of  outer  loops  from  32  to  10  or  15.  This  may  improve  the  convergence  of  the 
least  squares  error,  if  non-physical  solutions  are  also  rejected  at  each  iteration. 

(9)  Compare  the  background  subtracted  from  the  flux  with  the  latest  available 
natural  environment,  to  ascertain  the  correctness  of  the  background  subtrac- 
tion. If  significant  differences  are  found,  the  background  should  bo  re-added  to 
the  flux  and  the  natural  environment  subtracted. 

(10)  A  suggestion  proposed  by  Professor  D.G.  Shankland,  Air  Force  Institute  of 
Technology,  involves  a  significant  modification  to  the  fitting  method.  This 
method  involves  an  interpolating  function,  or  measure  of  "roughness"  of  the 
fit,  which  may  be  made  arbitrarily  smooth.  The  advantage  of  this  method  is 
that  no  limiting  assumptions  need  be  made  about  the  data.  The  chief  disad- 
vantage is  that  the  method  is  a  strictly  mathematical  method,  with  no  physical 
constraints.  Although  the  method  would  involve  extensive  repro^ramming,  it 
promises  worthwhile  information  about  the  flux  curves,  and  should  be  tried. 
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APPENDIX  A 

Demonstration  of  Zero  Diffusion  Current 

at  the  Magnetic  Equator  Tor  Half-Integer 

Ik'ssel  Function  Solution 

To  demonstrate  zero  diffusion  current,  one  must  show  that 

lim,_0Z>„  g'n{z)  =  0 

for  the  <7  =  0  case.  We  Know  that  from  (31): 


and  from  (16): 
and  from  (18): 


J-ir.(z)  =  y/2/(xz)   cos(~) 


Wc 


f-w  "  7^  ^  \te 


L)(*/2)(*/*( 


cos 


(<— >?f) 


/'(««) 


cos((2»+l)f  j). 


Hence,  the  derivative  of  the  eigenfunction: 

i'.W  -  - 


y/2jz~.s/ijl^(2n  +  1)) 


/'(*.] 


!n+l); 


sin((2n  +  1] 


2    2, 


Finally,  the  quantity  which  must  go  to  tero  is: 
D„g'n  =  - 


D„.  Vy7e  y/i/{7r*(2n+l))  ((2n  +  1)  Jr/(2ze)) 


■/'(«») 


sin 


(*♦.,£) 


The  term  in  square  brackets  is  a  finite  constant  for  any  n,  and  clearly  the  sine  goes  to  zero 
as  z  goes  to  tero.  Hence,  the  limit  is  zero. 


icn 


APPENDIX  B 


Derivation  of  Eigcnf unction  Derivatives  /^  and  /n 
To  find  the  derivative  of  the  steady  state  function,  f  x,  one  starts  with  (21): 


_  y  h.  I 


5V)3„(*W) 


If  the  Xn  are  assumed  (from  (19))  independent  of  z,  then  the  series  can  simply  be 
differentiated  term-by-term,  since  the  definite  integral  is  simply  a  constant  Tor  each  value 
of  n. 


dz 


""I        An         JO 


To  End  the  derivative  of  the  eigenfunction  j/„,  one  starts  with  (18): 


-=^a^/*4^^''/2',) 


Using  the  product  rule  for  differentiation: 


r.  =  sojo 


© ■'•(-©"■'■"'feo-ixri) 


'V{2-<r)/ze 
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APPENDIX  C 


The  Specific  Solution  of  the  Integral-Energy 

Omnidirectional  Flux  Equation  for  the  Half-Integer 

Bessel  Function  Solution 

The  integral-energy  omnidirectional  dux  Equation  (43),  which  is  the  energy-integrated 
form  of  Equation  (37),  is  represented  as  the  steady  state  term,  /oo,  plus  the  sum  of  exponen- 
tially decaying  eigenmodes  /o  through  /j: 

J>  e[v  =  VBo/b)  =  /oo  +  /o  +  A  +  h  +  h 


In  this  appendix,  each  term  is  shown,  as  derived  from  Equations  (37)  and  (43).  Note 
that  the  eigenmode  sum  should  be  infinite,  but  only  the  fundamental  mode  and  three  higher 
modes  are  shown. 

The  constants  in  each  term  are  defined  as  follows  (Ref.  39:55): 


T0  =  1  +  I  (in  2  +  y/3  )/(2>/3  H  »  1.3801730 

T,  =  ( Tv^/fl  ;=«  0.7404805 

K  =  (4/1  l)[r0  -  T,]  «  0.2328155 


The  "fitting  constants"  are  defined  as  in  Equation  (43)  or  in  Appendix  D: 

J  E 
a„  =  /     p-an(E)  dE 

JE 


=  23T.9  / 


J-VlSo/tf) 


(*/Ze)[x/Ze 


v/l-j/-(/?//i0) 


[To-iTo-W^y 


W*z.y 

»in[fr/2gc)([(l  -  y-)T0}/2  -  K{\  -  y"/'))]\ 

(sin  [(3r/2Ze)(((l  -  y2)T0)/2  -  K{\  -  y»/<))]) 

9 
(sin  [(5t/2Z«)(((1  -  r)T0)/2  -  K{\  -  y"/"1))]) 

25 
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(sin  [(7t/2,7c)(((1  -  y=)T„)/2  -  K{\  -  f  "/«))])  \ 

>  Ay 


49 


/o  =  -2™0'-"/n,)  / 


j-VSo/O 


v/sT^W  >/i-r(w//?o) 


To -(To-  n)t//< 


;in  (jt/(2Z.m(((1  -  y2)T0)/2  -  ff(l  -  y' ,/4))])  <iy 


/i  =  +2jra„er   f 


y/Uo/U 


y/2/Ze{Zx/Zc)\  \/l-y-(B/Bo)\\T0  -  (T0  -  T,  )yv' 


•s/Bt/B 


(sin  (3*/(2Zc))(((l  -  v-)T0)/2  -  A"(l  -  y'  «"))])  dy 

/2  =  -2Ta2e~,'/T»1  /  %/2/%(5*/Ze)  Vl  -  y-(AJ/^„)    Tn  -  (T,  -  T,)yV4  y 

(sin  (5jt/(2Z.))(((1  -  y?)T0)/2  -  A'(l  -  „' "«))])  ,<y 

/,  =  +2TO,e~(,/rj)  / 


■<Jb0/b 


y/l-y'2{B/Do)    To  -  (To  -  T,  )yv<  y 


(sin  [(7T/(2Ze))(((l  -  y2)T0)/2  -  /f  (l  -  y1  '/«))])  rfj, 

These  fitting  constants  arc  really  fuctions  of  energy  and  /./-shell. 

The  ratios  or  r0,  r( ,  r2,  and  r,i  should  be  (as  computed  from  the  ^igen  vn!u<'  ratio  of  the 
zeroes  of  the  function): 


r0  Rtf  9n 

ro  «a  25r? 
r0  £»  49  f3 


n  =(l/9)rn 
or  alternatively         ra  =  (l/25)n» 
r3  =  (l/49)r0 


APPENDIX  D 


The  General  Solution  of  the  Integral-Energy 

Omnidirectional  Flux  Equation  for  Bessel  Functions 

of  Arbitrary  Order,  as  a  Function  of  Sigma 

Substitution  of  (18),  (28),  and  {'29)  into  (27)  completes  the  general  expr^-.sion  for  the 
differential  energy  omnidirectional  flux,  J ,  with  Bessel  functions  of  order 

<7-  1 


—  a 


where  a  is  some  number  less  than  two: 


s/{2-o)lze 


+ 


(£> 


(l-<r)/2 


,-(t+1)/2 


L/l      "-t/fl 


,'t/j   K-vn 


l-l»/!| 


z.  V         2 


l-|»/3] 


-(<r/2) 


/)  4r)) 


-(»+l)/2        / 


-(<-/-! 


i-r(«/Bo)  r(y)  d{f) 

The  analogy  between  the  above  equation  and  Equation  (43)  is  seen  if: 
{Je   p2/00l(E)dE^    =  jE    p'S/D„t  dE 


and  if: 


\jE  p2an(E)      (c  -  t/(rn))L  =  j^  P-  dE  AnL{E) 

=  np-an(fi)c-t/ir"(n:))dE 


Hence,  (43)  is  simply  the  energy-integrated  form  of  (37). 
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APPENDIX  E 


Flow  Chart  for  Program  Electrofit 


Initial  guess  r0  (and  rt,  r*). 


I 


Numerically  integrate  functions  at 
each  data  point. 


Least  squares  (it  for  the  constants: 


Compare  fit  against  previous  "best' 
fit  using  least  squares  residual. 


better 


no  better 


Cheek  for  %  tries  with   no  im- 
provement 


Random  walk  r's 


no 


yes 


> 


— '    Optimize  yent 


Update  "best"  r's  and  a's. 


Cheek  for  7,1  better  fits 


yes 


no 


Print  results. 

Compute  fluxes. 

IMot   computed    flu-xcs  with   raw 

data. 
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